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SYMBOLS 

c  specific  heat  <BTU/lb  F> 

h  corrective  heat  transfer  coefficient  (BTU/hr  ft**2  F) 

hfg  heat  of  vaporization  (BTUVlb) 

rh  relative  humidity,  decimal 

t  time,  hours 

x  bed-depth  coordinate,  -ft 

G  dry  weight  flow  rate,  lb/hr  ft**2 

H  humidity  ratio,  lb/lb 

M  local  or  average  moisture  content,  dry  basis  (decimal) 

MR  moisture  ratio,  dimensionless 

P  pressure,  psia 

R  gas  constant,  ft  lb/lb 

S  cross-sectional  bed  area,  ft**2 

T  air  temperature,  F 

Tabs  absolute  temperature,  R 

V  velocity,  ft/hr 

P  dry  weight  density,  lb/ft**3 

9  product  temperature,  F 

i  bed  porosity,  decimal 
as  a  lower  case  symbol  beside  the  main  variable: 
air 


e 

equi 1 ibrium 

in 

inlet 

o 

at  time  t=o 

P 

product 

SYMBOLS CONT. 

saturated  vapor 
at  time  t 
uapor 
water 


INTRODUCTION 

During  the  past  -fine  years  at  Kansas  State  University,  research  has 
been  conducted  in  the  Mechanical  Engineering  Department  to  model  the 
performance  and  -feasibility  of  a  solar  powered  sorption  dehumidi  f  ier .  The 
dehumidifier  uses  a  rotating  bed  o-f  silica  gel.  This  desiccant  removes 
water  -from  the  air  stream  that  is  to  be  dehumidi-f ied.  Once  the  desiccant 
becomes  saturated,  it  is  regenerated  with  the  use  o-f  solar  heated  air.  As 
this  project  progressed,  a  practical  application  -for  the  use  of  the 
dehumidified  air  was  sought. 

The  sorption  dehumidi-f  ied  air  lends  itself  to  a  multitude  of 
applications  in  both  industry  and  agriculture.  Dehumidified  air  has  found 
use  in  the  production  of  foods,  environmental  countrol,  the  chemical  and 
pharmaceutical  industries,  and  the  drying  of  lumber.  Recently,  the  use  of 
solar  powered  dehumdified  air  to  dry  grains  has  become  a  topic  of  interest 
to  the  agricultural  community  concerned  wi th  the  fossil  fuel  depletion  and 
subsequent  rising  conventional  fuel  costs. 

Of  the  agricultural  crops  requiring  drying,  corn  uses  the  most  energy, 
is  the  largest  grain  crop  in  terms  of  total  production,  and  is  normally 
harvested  with  more  excess  moisture  than  any  other  grain  crop  (13).  It  has 
been  estimated  that  at  least  as  much  energy  is  used  in  drying  an  acre  of 
corn  as  is  used  for  all  the  other  farm  operations  necessary  to  grow  and 
harvest  that  acre,  including  operations  such  as  soil  preparation,  planting, 
cultivation,  and  harvesting  <16). 

The  present  research  applied  the  numerical  model  of  the  solar  energy 
powered  sorption  dehumidifier  to  the  drying  of  shelled  corn.  The 
objectives  of  this  research  were: 


1.  Decide  upon  a  type  of  grain  dryer  that  would  be  most 
applicable  in  ease  of  mathematical  modeling,  management, 
and  maintenance  for  research  purposes. 

2.  Once  a  certain  type  of  dryer  was  selected,  seek  out  the 
numerical  models  available  for  use. 

3.  Couple  the  selected  dryer  model  to  the  existing  sorption 
modei,  so  that  the  performance  of  the  complete  system 
can  be  predicted. 

4.  Uith  the  total  sorption  grain  drying  model  in  hand,  model 
<with  the  use  of  various  subroutines)  three  systems  for 
comparison  purposes: 

1.  Grain  dried  with  the  use  of  ambient  air, 

2.  Grain  dried  with  solar  heated  air. 

3.  Grain  dried  with  the  solar  sorption  dehumidified 


THE  DEHUMIDIFICATION  PROCESS 

Solid  sorption  can  be  accomplished  by  absorption  o-f  adsorption, 
depending  upon  whether  there  is  a  chemical  change  in  the  desiccant  during 
the  process.  In  the  absorption  of  water  vapor  on  a  desiccant,  the 
desiccant  undergoes  a  chemical  change  to  a  hydrate  state.  With  the 
addition  o-f  more  moisture,  absorption  will  cause  the  desiccant  to  dissolve 
into  solution.  In  adsorption,  the  desiccant  does  not  react  chemically  with 
the  condensed  water  vapor.  Sorbent  is  the  term  referring  to  the  desiccant, 
which  sorbs  moisture  -from  the  moist  air  stream  to  be  dehumidified.  The 
sorbate  is  the  substance  sorbed,  the  water  vapor  in  the  air. 

In  the  sorption  process  a  moisture  mass  transfer  takes  place  due  to  a 
water  vapor  partial  pressure  gradient  between  the  desiccant  and  the  air. 
This  process  will  continue  until  the  vapor  pressure  o-f  the  water  in  the 
desiccant  reaches  equi I ibrium  wi th  the  partial  pressure  of  the  water  vapor 
in  the  surrounding  air.  The  result  of  this  action  is  that  the  moisture 
content  of  the  air  decreases,  and  the  moisture  content  of  the  desiccant 
increases. 

Sorption  is  an  exothermic  process,  resulting  from  the  heat  of 
condensation  of  the  water  vapor  plus  the  heat  of  wetting.  This  last  term, 
the  heat  of  wetting,  occurrs  when  the  liquid  water  droplets  and  the 
desiccant  (silica  gel)  contact  one  another.  It  is  greatest  when  the 
desiccant  has  just  been  reactivated  and  tapers  off  as  the  sorbent  reaches 
saturation.  The  heat  of  condensation  and  the  heat  of  wetting  together  make 
up  the  heat  of  sorption. 

Solid  sorbents  generally  are  extremely  porous  "solid  foam",  with  large 
internal  surface  areas.   Silica  gel,  the  sorbent  used  in  this  study,  is 
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used  to  take  up  moisture  at  room  temperatures  from  an  air  stream  flowing 
through  it.  The  process  of  adsorption  by  solid  desiccants  is  reversible. 
Air,  when  warmed  by  passing  through  a  solar  collector,  has  an  increase  in 
its  moisture  carrying  capacity.  Figure  1.1  depicts  the  reason  for  this. 
At  a  lower  temperature,  Tl ,  the  maximum  amount  of  moisture  that  can  be 
"picked  up"  before  the  air  stream  reaches  saturation,  is  symbolized  byZ^Pl 
on  a  psychroraetr ic  chart.  However,  if  the  air  stream  is  sensibly  heated  to 
T2,  the  maximum  amount  of  moisture  that  can  be  held  by  the  air  stream,  is 
symbolized  by/JiP2.  ]n  practice,  the  raaximumO  P2  will  not  be  reached 
because  some  of  the  energy  in  the  air  is  lost  to  the  heating  of  the 
desiccant.  This  process  is  shown  as  a  dotted  line  on  Figure  1.1.  Because 
the  higher  air  stream  temperature  of  the  regeneration  process  causes  the 
air  to  be  further  from  saturation,  equilibrium  is  sought  between  the 
saturated  desiccant  and  the  air  stream,  and  water  vapor  is  transferred  from 
the  desiccant  to  the  air.  Silica  gel  is  especially  suited  to  this 
application  since  the  temperatures  needed  to  regenerate  this  desiccant  are 
those  attainable  by  flat  plate  collectors. 


Background 

The  solar  powered  sorption  dehumi di t i cat i on  system  in  this  study  was 
initially  developed  by  Singer  (28).  It  utilizes  a  desiccant  (silica  gel) 
to  dehumidify  air  through  adsorption,  with  solar  energy  to  regenerate  the 
desiccant.  Singer's  work  had  three  main  objectives: 

1.  To  select  an  optimum  solar  dehumi di f icat ion  technique  by  evaluat- 
ing the  numerous  dehumidif icat  ion  types  and  arrangements  with 
regard  to  their  application  to  the  use  o-f  solar  energy. 

2.  Develop  a  computer  model  simulating  this  solar  dehumidif icat i on 
process. 

3.  Construct  a  test  apparatus  of  this  system  so  an  evaluation  could 
be  made  o-f  the  computer  model. 

A  schematic  of  the  continuous  solar  air  dehumidif icat ion  system 
choosen  by  Singer  is  shown  in  Figure  2.1.  To  make  the  system  continuous, 
the  desiccant  wheel  is  rotated,  allowing  the  simultaneous  dehumi di t icat ion 
of  air  through  one  half  of  the  bed  and  desiccant  regeneration  in  the 
remaining  portion  of  the  bed.  The  flow  of  air  is  directed  through  the 
system  by  way  of  two  isolated  flow  paths:  the  process  flow  stream,  in 
which  air  will  be  dehumidified,  and  the  regeneration  flew  stream  which  is 
used  to  regenerate  the  desiccant. 

Silica  gel  was  selected  as  the  desiccant  material  to  be  used  in  the 
dehumidif ier .  The  reasons  for  this  were  threefold.  First,  silica  gel  is 
readily  available  commercially,  and  is  accepted  by  the  dehumidif icat ion 
industry.  Secondly,  there  is  publ i shed  mater i a]  available  regarding  the 
physical  characteristics  and  equilibrium  data  of  silica  gel  and  water  vapor 
(10)  (19)  (20).   Finally,  there  have  been  investigations  into  the  numerical 
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modeling  of   the   silica  gel  dehumi  d  i  f  i  cat  i  on  process  utilizing  -fixed  and 
rotating  beds  (24)  (24). 

Figure  2.2  shows  a  psychrometric  plot  o-f  the  -flow  streams  at  each 
state  during  the  dehumi  di-f  icat  i  ion  process.  Air  to  be  dehumidified 
(ambient  air)  enters  the  process  flow  stream  at  state  1  and  is 
adiabatical ly  dehumidi-f  ied  to  state  2.  The  dried  air  has  a  higher  dry  bulb 
temperature  resulting  from  the  heat  o-f  sorption  being  converted  to  sensible 
heat  in  the  process  airstream.  This  -flow  stream  then  passes  through  an 
optional  sensible  heat  exchanger  to  reduce  the  -final  dry  bulb  temperature 
to  state  3  i-f  desired. 

The  regeneration  flow  stream  begins  with  ambient  air  at  state  1.  The 
air  is  heated  to  state  4  by  passing  it  through  the  sensible  heat  exchanger, 
adding  the  heat  extracted  from  the  process  flow  stream.  Because  this  temp- 
erature is  too  low  for  desiccant  bed  regeneration,  the  flow  stream  is  then 
directed  through  the  flat  plate  solar  collector.  If  the  temperature  is 
still  not  sufficient  for  regeneration,  an  auxiliary  heater  energizes  to 
increase  the  temperature  to  state  5.  As  the  regeneration  air  stream  passes 
through  the  desiccant  bed,  the  bed  is  regenerated  and  the  air  stream  is 
humidified  to  state  6   where  it  is  released  to  the  atmosphere. 

The  model  adopted  by  Singer  to  simulate  the  rotating  desiccant  bed  was 
developed  by  Maclaine-Cross  and  Banks  (23).  Their  solution  involved  the 
transformation  of  mass  and  energy  conservation  equations  into  two  potential 
kinetic  wave  equations.  These  equations  were  linearized  and  solved 
numerically  to  provide  exit  air  stream  temperature  and  humidity.  Nelson 
(24)  applied  this  model  to  rotating  desiccant  wheels  and  published  a 
computer  program  to  apply  this  method  to  a  rotating  bed  of  silica  gel. 
Singer  modified  Nelson's  program  and  also  included  subroutines  to  model  the 
performance  of  the  solar  collector.   This  constituted  a  comprehensive 
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system  model  of  a  silica  gel  air  dehuriidif ier  with  solar  energy  powered 
regenerat  ion. 

A  physical  system  was  built  consisting  of  a  desiccant  wheel,  sensible 
heat  exchanger,  air  seals,  motor  drives,  dehumidif icat ion  housing  and 
ducting,  solar  collector,  blowers,  auxiliary  heaters,  and  instrumentation. 
Singer  performed  numerous  tests  to  compare  the  predicted  results  with  the 
experimental  results  obtained.  The  computer  program  was  found  to  be 
acceptable  in  predicting  the  performance  of  the  experimental  solar 
dehumi  di  f  i  er . 

Ananth  (2)  continued  the  verification  of  Singer's  work  by  collecting 
additional  data  from  the  experimental  apparatus.  He  studied  the 
performance  of  the  system  over  a  wide  range  of  parameters,  varying  air  flow 
rates,  temperatures,  humidity  ratios  and  peripheral  leakage  rates  for  the 
process  and  regeneration  streams.  Performance  curves  were  estabilished 
over  all  possible  ranges  for  the  variables  of  interest.  Ananth  also 
modified  the  computer  program  to  accept  varying  inlet  conditions  within  the 
time  frame  of  the  study. 

Following  Ananth's  work,  Atkinson  (4)  began  work  on  incorporating  TMY- 
SOLMET  weather  data  (3D  with  the  dehumi  di  f  i  er  model.  Since  the 
performance  of  a  solar  energy  powered  system  depends  on  weather  conditions, 
the  application  of  recorded,  "averaged"  weather  data,  would  allow  an  hourly 
estimation  of  the  performance  of  the  dehumi di f i cat i on  system.  The  TMY- 
SOLMET  weather  tape  provides  a  typical  meteorological  year  for  a  particular 
location.  The  typical  year  was  determined  using  statistical  methods  to 
select  a  typical  month,  for  each  of  the  twelve  calendar  months,  from  a  data 
bank  of  23  year's  observation.  The  tape  provides  weather  data  for  26 
locations  within  the  United  States.  Variables  such  as  extraterrestrial, 
beam,  diffuse,  and  engineering  corrected  radiation  Are     provided.   Other 
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variables  on   the   tape  are  dry  bulb  temperature,  dew  point  temperature, 
atmospheric  pressure,  and  absolute  humidity. 

Once  the  performance  of  the  dehumi  di  f  i  cat  ion  process  was  modified  to 
accept  location  weather  data,  Atkinson  began  work  on  an  economic 
■feasibility  study.  The  method  choosen  was  the  PI,  P2  method  outlined  by 
Duffie  and  Beckman  <12>.  This  method  relates  present  worth  factors  for  the 
life  cycle  costs,  which  is  merely  the  present  worth  of  all  costs.  This 
method  is  applied  to  both  the  conventional  and  solar  system  to  provide  a 
foundation  for  comparision.  The  reader  is  encouraged  to  consult  Atkinson's 
thesis  (4>  for  further  details. 

Throughout  the  course  of  studying  the  sorption  dehumidif ier,  the 
investigators  Singer,  Ananth,  and  Atkinson  modified  the  computer  model  many 
times,  for  various  reasons.  It  would  be  appropriate  at  this  time  to 
summarize  the  state  of  this  program  and  it's  subroutines  when  the  present 
investigator  began  work.  The  main  program  served  as  the  basis  for  reading 
the  location  weather  data  and  the  calculation  of  the  economic  analysis. 
There  were  problems  with  the  reading  of  the  weather  data,  because  an 
algorithm  had  not  been  developed  which  would  allow  the  crossing  of  month 
boundaries  in  reading  the  tape.  The  first  day  of  the  year  was  also 
unattainable.  The  main  program  allowed  for  the  input  of  the  economic 
parameters  and  functioned  as  the  mechanism  for  calling  the  subroutines. 
There  were  five  subroutines  associated  with  the  dehumidifier  model: 
SYSTEM,  COLPER,  HUMID,  ALFAV,  EGU. 
SYSTEM 

The  subroutine  SYSTEM  read  data  for  the  following  parameters: 
1.   The  program  option,  IFLAG,  indicated  whether  the 

temperture  of  the  regeneration  air  into  the  dehumid- 
ifier was  set  to  a  minimum  temperature  reached  by 
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the  use  of  an  auxiliary  heater  <H2),  or  if  the 
entering  air  temperature  was  to  be  set  to  the 
outlet  temperature  o-f  the  solar  collector  (#1). 

2.  Collector  data  such  as  the  air  volume  flow  rate  (1/s), 
optimum  slope  of  the  collector  <see  appendix  A), 
latitude  and  longitude  were  read. 

3.  Dehumidif  ier  data  including  duct  temperature  drops,  the 
geometrical  properties  o-f  the  silica  bed,  the  revolutions 
per  second  o-f  the  desiccant  wheel,  the  minimum  regeneration 
temperature  (i-f  using  option  *2,  IFLAG),  leakage  rates, 

and  heat  exchanger  efficiency  were  also  provided. 

SYSTEM  called  the  subroutines  COLPER  and  HUMID.  Output  from  this 
subroutine  was  provided  in  the  form  of  an  echo  printing  of  the  dehumidif ier 
operating  parameters.  An  hourly  temperature  and  humidity  map  of  the 
process  side  of  the  dehumidifer,  and  an  hourly  temperature  and  humidity  map 
of  the  regeneration  side  was  also  provided.  A  total  heat  balance  map  of 
the  dehumidif ier,  heat  exchanger,  ducting  to  and  from  the  collector,  the 
collector,  the  auxiliary  heater,  and  the  surplus  solar  heat  available  could 
also  be  found  in  this  subroutine. 
COLPER 

The  subroutine  COLPER  mode  1 ed  the  performance  of  the  solar  collector. 
It  was  originally  developed  by  Singer  and  based  on  equations  in  Duffie  and 
Beckman  (12).  The  subroutine  determined  the  useful  heat  gain  in  the 
collector  by  calculating  the  amount  of  radiation  striking  the  collector. 
It  did  this  by  utilizing  ambient  conditions,  collector  orientation,  time  of 
day,  and  time  of  year.  Knowing  the  incident  radiation  on  the  surface  of 
the  collector,  the  subroutine  proceeded  to  obtain  the  collector  efficiency 
given  ambient   conditions,   inlet   collector   conditions,   and   loss 


13 

coefficients.   Once  the  efficiency  was  calculated,  the  collector  useful 
heat  gain  and  outlet  temperature  were  determined. 

The  COLPER  subroutine  also  established  the  amount  of  heat  lost  through 
the  ductwork  from  the  collector  to  the  dehumidi f ier,  given  estimated 
temperature  drops.  If  the  temperature  returning  to  the  dehurai di-f i er  was 
below  the  required  inlet  regeneration  temperature,  the  COLPER  subroutine 
found  the  amount  of  auxiliary  heat  required  to  assure  proper  regeneration 
temperature.  The  amount  of  heat  available  for  storage  was  also  determined. 
The  collector  modeled  in  this  section  utilized  the  manufacturer's 
performance  curves  for  a  Solaron  Series  3000  collector,  but  could  be 
modified  for  other  manufactured  collectors. 
HUMID 

THE  HUMID  subroutine  analytically  modeled  the  simultaneous  heat  and 
mass  transfer  occuring  within  the  desiccant  bed.  The  model  developed  by 
MacLame-Cross  and  Banks  (23.)  and  computerized  by  Nelson  (24),  predicted 
the  performance  of  a  rotating,  silica  gel,  desiccant  bed.  The  subroutine 
provided  process  and  regeneration  outlet  temperatures  and  humidity  ratios. 
Please  refer  to  Singer  (28)  and  Ananth  (2)  for  further  details. 
ALFAV 

ALFAV  was  a  support  subroutine  for  the  HUMID  subprogram.   It  supplied 
the  equilibrium  properties  of  the  moist  silica  gel  and  the  air  water  vapor 
mixture.  Nelson's  Thesi s  (24)  provides  an  indepth  explanation  of  the  exact 
function  of  this  subroutine,  please  consult  it  for  further  details. 
E6V 

EGV  (Equilibrium  Gamma  Yalues)  was  also  a  support  subroutine  for  the 
HUMID  subprogram.  It  determined  the  incremental  steps  that  would  be  taken 
in  determining  the  temperature  and  humidity  ratios  during  the  computational 
procedure.   Consult  Nelson   (24)  or  Atkinson  <4)  for  clarification  of  the 
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purpose  of  this  segment. 

All  o-f  the  subroutines  and  the  main  program  had  to  be  able  to  interact 
with  variables  of  common  interest.  This  was  done  through  the  use  o-f  large 
common  blocks  that  allowed  the  trans-fer  of  information  from  one  subprogram 
to  another.  Examples  included  wi thin  the  common  blocks  were:  the  weather 
data,  the  inlet  and  outlet  conditions  of  the  major  components  of  the 
system,  counters,  and  properties  of  state,  location,  and  geometry. 

This  thesis  completes  the  next  step  in  the  development  of  the 
application  of  this  sorption  dehumidifier  model.  The  application  is 
directed  toward  grain  drying. 
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LOAD  MODEL 

The  amount  of  moisture  in  grain  has  an  effect  on  its  performance  for 
such  processes  as  harvesting,  storage,  and  germinating.  If  the  grain  is 
too  "wet"  it  will  provide  an  environment  in  which  molds  and  insects  will 
thrive.  If  the  grain  is  "over-dried",  it's  ability  to  germinate  can  be 
adversely  effected.  In  practice,  grain  will  not  be  harvested  at  moisture 
contents  greater  than  35%  and  grain  drying  will  not  usually  continue  for 
moisture  contents  below  10%  (wet  basis)  (21). 

The  grain  moisture  contents  in  Table  3. A  are  those  recommended  for  the 
safe  storage  of  grain.  The  length  of  time  that  crops  can  be  stored  varies 
with  the  moisture  content  and  type  of  crop.  To  store  a  crop  for  5  years 
the  moisture  content  should  be  approximately  2M  below  the  moisture  content 
that  is  considered  safe  for  a  1  yar   storage  (8). 


Tabic  3. A 
Moisture  Content  During  Harvest  and  for  Safe  Storage,  Percent,  w.b. 


Requ  i  red 

for 

Safe 

Maximum 

Opt  imum 

at 

Storage 

During 

Harvest 

for 

Usual  when 

Cereal 

Harvest 

Min  imum 

Loss 

Harvested 

for  1  yr 

fi 

ar  5  yr 

Barley 

30 

18-20 

10-18 

13 

11 

Corn 

35 

28-32 

14-30 

13 

10-11 

Oats 

32 

15-20 

10-18 

14 

11 

Rice 

30 

25-27 

14-25 

12-14 

10-12 

Rye 

25 

\6-20 

12-18 

13 

11 

Sorghum 

35 

30-35 

10-20 

12-13 

10-11 

Wheat 

38 

18-20 

9-17 

13-14 

11-12 

From  Drying  Cereal  Grains,  BrooKer,  Bakker-Arkema,  and  Hall,  1982. 
Sources:   C.  U.  Hall  (1957>j  D.  W.  Hall  (1970)j  Matz  (1969);  Sinha  (1973). 


The  objective  in  grain  drying  is  to  reduce   the  moisture   content   so 
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that  spoilage  will  not  occur  be-fore  the  grain  can  be  used.  Drying  a  high- 
moisture  grain  at  an  original  moisture  content,  (Mo) ,  to  a  -final  moisture 
content,  <M-f),  can  be  carried  out  ouer  a  long  period  o-f  time  if  a  low 
drying  air  temperature  is  used;  less  time  is  required  when  a  higher  air 
temperature  is  used. 

Temperature  is  not  the  only  parameter  that  influences  the  time 
required  to  reduce  the  moisture  content  o-f  a  grain.  Air  relative  humidity, 
airflow,  initial  moisture  content  and  -final  moisture  content  all  dictate 
the  amount  o-f  time  that  will  be  required  to  reach  the  desired  condition  o-f 
the  grain.  The  usual  ranges  of  airflow  and  temperature  are  listed  in  Table 
3.B  (8). 


Tabic  3.B 
Usual  Range  of   Air-flows  for  Drying  Systems,  CFM/BU 


Aeration  1/50  -  l 

Natural  Air  (unheated)  2-5 

Layer  Drying  2-10 

Heated  Air  < 130  -  500  F)  30  -  100 


A  typical  drying  curue  is  shown  in  Figure  3.1.  The  reader  will  notice 
the  characteristic  rapid  rate  o-f  drying  initially,  and  the  slower  rate,  as 
the  total  time  o-f  the  drying  increases  (21).  As  the  moisture  content  gets 
closer  to  equilibrium  conditions  with  the  air,  the  drying  becomes  so  slow 
that  is  can  hardly  be  detected.  This  description  applies  to  the  drying  of 
grain  of  any  kind  that  is  fully  exposed  to  an  atmosphere  of  constant 
temperature  and  humidity.  Generally,  the  rate  o-f  drying  is  faster  if  the 
initial  moisture  content  was  higher.  High  drying  temperatures  and  low  air 
absolute  humidity  will  also  result  in  -faster  drying  rates. 
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Figure  3.1 
Typical  Drying  Curve 
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At  this  point,  some  important  concepts  in  the  understanding  of  grain 
drying  will  be  de-fined;  Equi  1  ibrium  Moisture  Content  <EMC>,  Moisture  Ratio 
<MR>,  and  the  reporting  of  moisture  content  on  a  wet  or  dry  basis  <w.b.  and 
d.b.,  repect iuely) . 
Equilibrium  Moisture  Content 

The  equilibrium  moisture  content  determines  the  minimum  moisture 
content  to  which  grain  can  be  dried  under  a  given  set  of  drying  conditions. 
It  is  defined  as  the  moisture  content  of  the  material  after  it  has  been 
exposed  to  a  particular  environment  for  an  infinitely  long  period  of  time. 
EMC  depends  upon  the  humidity  and  temperature  conditions  of  the 
surroundings  as  well  as  the  maturity  and  type  of  grain  (8). 
Moisture  Ratio 

The  moisture  ratio  is  a  quantity  often  referred  to  in  grain  drying 
nomenclature.   It  is  defined  as: 

MR  =  <  M  -  Me)  /    (Mo  -  Me; 

where  - 

MR  =  moisture  ratio. 
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M  =  current  moisture  content  of  the  grain. 

Me  =  desired  moisture  content  of  the  grain  at  the  end 

o+  drying. 

Mo  =  moisture  of  the  grain  at  the  start  of  drying. 

Moisture  Content  on  a  Met  Basis 

The  moisture  content  o-f  a  grain  is  used  as  a  measure  o-f  maturity  and 

quality.   The  elevators,  or  market  place,  are  usually  interested  in  the 

moisture  content  on  a  wet  basis.  This  is  because  the  denominator,  the  wet 

weight,   is  the  quantity  obtained  when  the  truck  loaded  with  the  grain  is 

weighed.  The  definition  of  a  wet  basis  moisture  content  is  <8): 

Mw  wet  basis  =  u z d  (100) 

w 
where:  w  =  wet  weight    d  =  dry  weight 

Mw  =  moisture  content  on  a  wet,  percent  basis 

Moisture  Content  on  a  Dry  Basis 

The  moisture  content  of  a  grain  on  a  dry  basis  is  used  in  many  engineering 

calculations.  Definition: 

Md  dry  basis  =  u = d  (100) 

d 
where:  Md  =  moisture  content  on  dry  basis,  percentage 

For  example,  if  we  are  told  that  grain  has  a  moisture  content  of  25X, 
we  might  expect  this  to  mean  that  100  lb.  of  grain  contains  25  lb.  of 
waterj  this  is  correct  if  the  moisture  content  is  expressed  in  per  cent  wet 
basis.  On  the  other  hand,  it  is  just  as  reasonable  to  assume  that  the  25'A 
moisture  content  is  expressed  on  a  dry  basis;  in  this  case  100  lb.  of  grain 
contains  20  lb.  of  water  and  80  lb.  of  dry  matter  since  20  is  25X  of  30. 
Accordingly,  when  the  moisture  content  is  reported  as  percentage  it  is 
necessary  to  have  an  understanding  as  to  which  basis  is  used.  One  is 
neither  more  correct  nor  more  logical  than  the  other  (21). 

In  selecting  a  grain  drying  mathematical  model  to  be  coupled  with  the 
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solar  sorption  dehumidif ier  model,  a  dryer  system  had  to  -first  be  selected. 
There  are  several  types  of  dryer  designs  that  were  considered.  They 
basically  fell  into  two  broad  categories:  those  that  dry  grain  in  batches, 
and  those  that  dry  grain  as  it  -flows  continuously  through  the  equipment. 


Batch  Drying 

Fixed  Bed  Drying 

This  type  of  system  can  use  a  bin  with  diameter  -from  a  few  -feet,  to 
many  feet.  It  is  usually  deep  (up  to  Id  feet)  and  has  a  perforated  -floor 
through  which  a  relatively  low  air-flow  rate  is  provided.  The  inlet  air  can 
be  ambient  or  heated.  The  grain  can  just  "set  there"  or  it  can  be  stirred. 
Characteristic  of  fixed  bed  drying  is  the  phenomenon  of  "drying  zones". 
Figure  3.2  is  a  schematic  of  a  fixed  bed  dryer.  As  drying  air  moves 
upward,  there  is  an  exchange  of  moisture,  from  grain  to  air,  in  a  finite 
depth  or  zone  of  grain.  At  the  start  of  the  drying  process  the  zone  is 
located  at  the  bottom  of  the  dryer.  As  the  drying  continues  the  zone  moves 
upward,  and  the  grain  has  been  dried  to  EMC  when  the  zone  reaches  the  top 
of  the  bed. 

There  are  several  different  types  of  batch  systems  that  will  not  be 
discussed  further  because  they  are  too  difficult  to  apply  as  a  mathematical 
model  and  they  require  careful  management  in  use.  Brooker,  et  a).,  (8) 
provide  a  thorough  discussion  on  the  subject  in  the  chapter  titled  "Grain 
Drying  Systems" . 
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Continuous  Flow  Dryers 

Cross-flow  Dryers 

Figure  3.3  is  a  sketch  of  a  typical  cross-flow  dryer,  in  which  grain 
flows  from  a  wet-grain  holding  bin  at  the  top,  down  the  columns,  and  is 
discharged  at  the  bottom.  The  upper  portions  o-f  the  columns  are  drying 
sections  and  the  lower  portions  are  cooling  sections.  A  metering  device 
and  temperature  sensor  are  used  to  regulate  the  flow  rate  of  the  grain. 
The  name  "Crossflow"  comes  from  the  fact  that  the  flow  of  heated  air  is 
perpendicular  to  the  flow  of  the  grain  (8). 
Concurrent  Dryers 

In  concurrent  dryers  the  air  flows  in  the  same  direction  as  the  grain. 
Figure  3.4  is  a  schematic  of  such  a  dryer.  There  are  no  large  central 
plenums  used  in  this  dryer  as  are  used  in  the  crossflow  dryer.  The 
characteristics  of  a  concurrent  dryer  are  a  series  of  small  ducts  used  to 
introduce  the  air  and  to  also  provide  for  the  escape  of  the  air.  Heated 
air  is  forced  into  the  upper  row  of  ducts,  and  cool  air  is  forced  into  the 
lower  ducts.  Both  the  heated  and  cool  air  exhaust  through  a  duct  that 
extends  across  the  dryer.  Ulet  grain  is  preheated  as  it  moves  downward 
along  the  ducts  carrying  heated  air  with  it.  After  the  grain  passes  the 
lower  edge  of  the  ducts,  the  air  and  the  grain  move  in  the  same  direction. 
The  hottest  air  enters  the  wettest  grain  and  is  quickly  cooled.  The  grain 
temperature  is  lower  than  the  air  temperature  at  this  point  in  the  flow, 
because  a  high  rate  of  evaporation  is  taking  place.  As  the  grain  continues 
to  move  downward,  its  temperature  increases  and  then  decreases  along  with 
that  of  the  drying  air  (3). 
Counterflow  Dryers 

In  counterflow  dryers  the  air  flows  in  the  opposite  direction  of  the 
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movement  of  the  grain.  This  method  can  be  visualized  by  considering  a  bin 
which  holds  a  group  of  thin  layers  of  grain.  The  drying  air  blows  upward 
through  the  layers.  At  each  time  interval,  a  new  layer  is  placed  on  the 
top  of  the  stack  and  a  layer  removed  from  the  bottom.  See  Figure  3.5  for 
an  example  counterflow  dryer  (28)  <8) . 

One  of  the  biggest  drawbacks  of  a  continuous  flow  dryer  for  research 
purposes  is  the  extra  equipment  needed  for  material  handling,  and  the 
problems  that  could  arise  in  the  maintenance  and  management  of  the 
necessary  moving  mechanical  parts.  Usually,  there  is  assoc iated  wi th  the 
continuous  flow  systems  a  considerable  amount  of  expensive  handling 
equipment  (augers,  metering  devices,  temperature  sensors).  In  fact,  the 
capacity  of  such  a  system  may  actually  be  limited  by  the  capacity  of  the 
auger  system.  Maintenance  of  all  mechanical  equipment  would  also  have  to 
be  considered. 

Another  problem  that  arises  when  considering  the  use  of  a  continuous 
flow  system  is  the  higher  temperatures  needed  (200-500F).  Since  there  is 
less  contact  time  during  the  process  between  the  drying  air  and  the  grain, 
the  higher  air  temperatures  are  needed  to  obtain  the  same  results  as  in  a 
batch  dryer.  This  results  in  a  inefficent  use  of  the  energy  in  the  hot  air 
because  a  good  portion  of  that  energy  is  used  in  heating  the  grain.  The 
application  of  a  solar  collector  to  heat  the  air  entering  a  continuous  flow 
system  would  present  several  problems.  Costs  of  collector  systems  to 
provide  high  temperatures  are  considerably  greater  than  for  lower 
temperature  systems.  Collection  efficiencies  are  reduced  in  high- 
temperature  collectors  unless  expensive  measures  are  taken  to  limit  heat 
losses  (13). 

On  the  other  hand,  a  batch-type,  fixed  bed  dryer  has  the  following 
advantages: 
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1.  The  grain  can  be  harvested  at  any  rate  desired. 

2.  Simple  management. 

3.  Minimum  grain  handling. 

4.  The  energy  in  the  drying  air  is  eff  ic  iently  used. 

5.  The  grain  is  not  ovei — dried. 

6.  The  low-temperature  air  causes  a  minimum  of  cracking. 

7.  The  temperatures  needed  -for  batch  drying  are  obtainable  with 
•flat  plate  collectors. 

Although  a  batch  dryer  may  not  dry  as  quickly  as  a  continuous  flow  dryer, 
it  can  provide  an  excellent  basis  for  comparisons  of  feasibility  and  performance 
when  varying  the  drying  air  temperature  and  humidity.   Therefore,  the  dryer 
choosen  for  this  study  was  a  fixbed  bed,  batch  dryer.   Once  the  dryer  type  was 
selected,  a  mathematical  model  for  the  fixed  bed  dryer  was  investigated. 
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Review  of  Grain  Drying  Models 
Research  was  done  to  determine  what  types  of  drying  models  were 
available  -for  use.  Hukill  (21)  in  1954  proposed  a  model  that  expressed 
analytically  the  moisture  o-f  the  grain  as  a  function  o-f  bed  height  and 
time.  More  recently,  Thompson,  et  al ,  (27)  developed  mathematical  grain 
drying  models  for  concurrent  flow,  crossflow,  and  counterflow  grain  dryers. 
These  models  were  empirical  in  nature  and  apply  only  to  corn  drying.  Grain 
drying  models  based  on  laws  of  heat  and  mass  transfer  lead  to  complicated 
systems  of  equations  that  can  only  be  solved  with  the  use  o-f  large 
computers  (but  can  more  generally  be  applied  to  other  biological  products). 
Bloome  and  Shove  (7),  and  Barre,  et  al ,  (4),  made  a  number  of  assumptions 
in  their  analyses  to  simplify  the  solution  of  the  the  heat  and  mass 
transfer  drying  equations.  Their  assumptions  have  been  claimed  only 
partially  valid  by  colleagues  (8).  The  most  general  application  of  the 
fundamental  laws  of  heat  and  mass  transfer  to  drying  biological  products 
was  made  by  Bakker-Arkema,  et  al ,  (5). 

The  grain  drying  model  developed  by  Bakker-Arkema,  et  al ,  at  Michigan 
State  University  was  the  model  choosen  in  this  study.  The  reasons  for  this 
were: 

1.  The  equations  of  mass  and  heat  transfer  are  general 
in  nature,  and  therefore  could  be  applied  to  other 
biological  products  when  an  appropriate  thin-layer 
equation  could  be  determined. 

2.  Fully  documented  computer  listings  for  the  modeling 
of  fixed  bed,  crossflow,  concurrent,  and  counterflow 
driers  are  included  (5).   If  the  type  o-f  dryer  was 
later  changed  in  some  further  study,  the  models  would 
be  available. 


28 

3.  The  grain  drying  models  could  be  considered  a  black 
box  which  only  required  certain  input  in -format  ion  in 
order  to  -furnish  the  desired  output. 

4.  The  models  developed  at  MSU  are  capable  o-f  predicting 
the  performance  o-f  -fixed  bed,  crossflow,  concurrent, 
and  counter-flow  driers  to  within  10X  of  experimental 
drying  rates  and  temperatures  (25). 

The  -fixed-bed  model  is  based  on  ideas  o-f  Schumann  (1929),  Van  Arsdel 
(1955),  Klapp  <19<S1),  and  Bakker-Arkema,  et  al  (1967).  The  highlights  of 
the  model  as  described  in  Drying  Cereal  Grains  (8)  will  be  detailed  in  the 
following  pages. 

Initially  an  elemental  bed  volume  is  drawn  as  shown  in  Figure  4.1. 
Energy  and  mass  balances  are  written  on  a  differential  volume  (Sdx).  The 
four  unknowns  in  this  system  of  equations  »re:  M,  the  average  grain  kernel 
moisture  content;  H,  the  humidity  ratio  of  the  air;  T,  the  air  temperature; 
and  6,  the  kernel  temperature.  Four  equations  must  be  derived  to  solve  for 
the  four  unknowns. 

The  four  equations  required  to  solve  this  model  follow: 

I.)  Enthalpy  of  the  Air 

Energy  out  =  energy  in  -  energy  transferred  by  convection 

(Gaca  t  GacvH>«<T  ♦OT/  3x)dx)Sdt=(Gaca  +  GacvH)STdt  -  haSdx(T-8)dt 

2.)  Enthalpy  of  the  Product 

energy  transferred  =  change  in  internal  product  energy  -  energy  for  evaporation 

haSdx(T-  e>dt  = 

(Ppcp   +    PpcwH>Sdx(3  9/ 3t)dt 

-Chfg   +   cv<T   -   &)]Ga('3H/'3x)dxSdt 
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3.)      Humidi  ty  of    the  Air 
moisture    transferred  ■  moisture    in   -  moisture   out 

CpSdx(3M/3t)dt  =  GaSHdt  -  GaS(H  ♦<  3  H/3x)dx)dt 

t 

4.)  Moisture  Content 

<9M/"9t>  «  an  appropriate  thin  layer  equation. 

These  equations  constitute  the  simulation  model  for  the  fixed  bed 
grain  dryer.  Since  an  analytical  solution  to  the  system  of  equations  is 
impossible,  the  differential  equations  are  solved  by  finite  difference 
techniques. 

The  initial  and  boundary  conditions  of  the  corn  and  the  drying  air 
must  be  known  in  order  to  solve  the  equations.  The  known  values  must 
i  nc  1  ude : 

1.  The  initial  temperature  and  moisture  content  of  the  grain. 

2.  The  initial  temperature  and  humidity  of  the  drying  air. 
The  specific  boundary  conditions  for  the  fixed  bed  dryer  are: 

T<o,t>  ■  T(inlet) 
*<x ,o)  =  9( initial) 
H<o,t)  ■  H(inlet) 
M<x,o)  -  M(initial) 
Where  T  =  temperature  of  the  air,  S  =  temperature  of  the  grain,  H  = 
humidity  ratio  of  the  air,  M  ■  moisture  content  of  the  grain. 

The  program  uses  the  following  sequence  in  solving  the  differential 
equations: 

1.)  input  data 

2.)  initialize  arrays 

3.)  evaluate  constants  used 

4.)  solve  the  differential  equations 

5.)  ouput  when  appropriate 
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To  solve  the  model   equations,  the  values  of  T,  H,  9,  and  M  must  be 

speci-fied  at  each  position  within  the  bed  before  the  dryer   is  started. 

Bakker-Arkema,   et  al ,  found  the  following  method  to  be  the  most  stable  and 

re  I i  able: 

T(x,o)  =  T(inlet) 
6<x,o>  ■  e< initial  i 
H(x,o)  =  H(inlet) 
M<x,o>  =  hKini  tial) 
The  conditions  above  are  physically  incorrect  because  they  assume   the 
■first  blast  of  drying  air  displaces  all  air  within  the  dryer  without  heat 
or  mass  transfer.  The  bed  is  initialized  this  way  strictly  for  stability. 
Also  covered   in  this  initialization  process  is  the  grain  temperature  next 
to  the  air  inlet.  This  temperature  is  set  equal   to  the  average  of   the 
inlet  air  and  the  initial  grain  temperature: 
«<o,o>  ■  (T<inlet>  *   e< ini t i al >>/2 
The   third  initialization  step   is  to  solve  for  all  values  of  the 
absolute  humidity,  moisture  content,  air  temperature,  and  grain  temperature 
for  each  position  in  the  bed  for  the  first  time  through   the  calculations. 
Figure  4.2  shows  the   indexing  scheme  for  the  pertinent  values  in  this 
study. 

The  portion  of  the  computer  program  that  models  the  fixed  bed  dryer  is 
made  up  of  a  main  program  FIXBED  and  three  subroutines;  LAYEQ2,  READYT, 
ZEROIN.  There  are  also  functions  called  upon  in  the  calculations  such  as 
EMC,  and  all  the  functions  included  in  the  SYCfWRT  PACKAGE.  Common 
properties  such  as  the  heat  capacity  of  air  and  water,  atmospheric 
pressure,  and  the  bulk  density  of  the  grain  were  made  available  to  all  the 
subroutines  and  functions  through  the  use  of  a  BLOCKDATA  FORTRAN  option. 
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FIXBED 

This  is  the  main  program  for  the  numerical  modeling  of  the  -fixed  bed 
grain  dryer.  It  is  within  this  segment  that  all  the  initialization  takes 
place  and  the  actual  depth  calculations  occur.  All  constants  needed  are 
calculated,  and  ail  output  -for  the  dryer  originates  here.  It  also  serves 
as  the  basis  -for  calling  needed  subroutines  and  functions. 

Function  EMC 

EMC  is  a  function  subroutine  that  computes  the  Equi I i brium  Moisture 
Content  o-f  corn  -from  a  given  relative  humidity  and  temperature.  The 
equations  used  sr^: 

For  temperatures  less  than  235F  the  DeBoer  equation  -(see  appendix  B) . 
For  temperatures  greater  than  235F  the  Thompson  equation  -(see  appendix  C) . 

LAYEQ2 

The  subroutine  LAYEQ2  is  the  thin  layer  drying  equation  that 
calculates  the  moisture  content  o-f  the  grain  as  it  varies  with  time.  There 
are  thin  layer  equations  that  apply  to  different  grain  temperature  ranges. 
The  equations  used  in  this  study  were  the  Subbah  equation  (see  appendix  D) 
■for  grain  temperatures  below  80F,  and  the  Troeger  equation  (see  appendix  E> 
for  grain  temperatures  between  80F  and  140F.  The  question  arose;  could  the 
Troeger  equation  be  used  as  a  less  accurate,  but  relatively  reliable 
predictor  o-f  moisture  contents  -for  temperatures  less  than  80F?  The 
benefits  of  using  one  thin  layer  equation  would  allow  the  temperature 
boundary  of  80F  to  be  crossed  in  a  study  time  frame.  This  would  eliminate 
the  need  for  complicated  checks  within  the  bed  to  determine  which  equation 
to  use.  If  the  Troeger  equation  was  found  to  be  unsatisfactory  in 
predicting  moisture  contents  for  grain  temperatures  below  80F,  how  was  the 
situation  to  be  handled  if  both  equations  were  required  at  different 
posit  ions  within  the  bed?  To  answer  this  question,  a  trial  run  was  made. 
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Ambient  air  for  Omaha,  NE  was  used  as  the  dryer  inlet  air  -for  the  day  of 
October  20.  The  simulation  had  the  air  enter  a  fixed  bed  dryer  of  five 
feet  in  height  and  one  square  foot  in  cross  sectional  area.  The  CFM  was 
set  at  120,  and  the  grain  temperature  was  initially  set  at  40F.  The  initial 
corn  moisture  content  was  33ft  dry  basis.  There  were  15  increments  per  foot 
in  the  calculations.  In  the  first  simulation  the  Troeger  equation  was 
used;  in  the  second  the  Subbah.  The  results  can  be  found  in  Table  4. A  and 
Graph  HI.  The  Troeger  equation  was  found  to  be  unacceptable  in  predicting 
moisture  content  in  grain  below  80F.  Therefore  ,  if  the  study  crossed  the 
boundary  of  30F,  different  thin  layer  equations  would  have  to  be  employed. 

READYT 

This  subroutine  is  a  support  program  for  LAYEQ2  ,  it  makes  preliminary 
checks  and  calculations  for  the  thin  layer  equation  and  calculates  the 
moisture  ratio. 

SYCWtRT  PACKAGE 

The  SYCHART  group  of  19  function  subprograms  and  one  subroutine  are  a 
numerical  modeling  of  the  psychrometr ic  chart.  Lerew  (22)  programmed  this 
collection  of  theoretical  and  empirical  psychrometr ic  equations  as  a  set  of 
interconnected  FORTRAN  subprograms.  Dry-bulb,  wet-bulb  and  dew-point 
temperature,  humidity  ratio,  relative  humidity,  vapor  pressure,  enthalpy 
and  specific  volume  equations  are  included  in  the  model.  If  any  two 
independent  properties  of  moist  air  are  known,  it  is  possible  to  find  the 
remaining  properties. 

2ER0IN 

2ER0IN  is  a  root  finding  technique.  It  is  based  upon  an  algorithm 
which  is  a  mixture  of  linear  interpolation,  extrapolation  and  bisection. 
ZEROIN  will  issue  a  warning  when  there  are  no  roots  or  multiple  roots 
between  the  initial  guesses.   If  this  error  message  appears,  it  normally 
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TABLE  4. A 


VALUES  OF  MOISTURE  CONTENT 
SIMULATED  USING  THE  TROEGER  AND  SUBBAH  EQUATIONS 


input  conditions:  product  temperature  60F 
moisture  content  33X 
120  CFM 
ambient  air  is  the  drying  air 


MOISTURE  CONTENTS 
AT  THE  SECOND  FOOT  IN  THE  BED 
HOUR  TROEGER  SUBBAH 


I  .3269  .3333 
3  .3269  .3333 
5  .3270  .3333 
7  .3283  .3340 
'  .3305  .3347 

II  .3292  .3347 
13  .3254  .3347 
15  .3173  .3347 
17  .3080  .3347 
1?  .3019  .3347 
21  .3003  .3347 
23  .2996  .3347 
25  .2985  .3347 
27  .2977  .3347 
29  .2970  .3347 
31  .2958  .3347 
33  .2932  .3347 
35  .2862  .3347 
37  .2755  .3347 
39  .2660  .3347 
*»!  .2591  .3347 
43  .2545  .3347 
45  .2502  .3347 
47  .2468  .3347 
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indicates  instability  requiring  a  larger  value  -for  the  number  of  layers 
calculated  per  -foot. 

MODIFICATIONS  TO  THE  FIXBED  PROGRAM 

In  order  to  apply  this  model  to  the  existing  solar  sorption 
dehumidif ier ,  the  program  had  to  be  modified  to  accept  varying  inlet 
conditions  of  air  absolute  humidity  and  temperature.  Problems  also  arose 
in  the  repetition  of  some  variable  names  used  in  both  the  fixed  bed  program 
and  the  existing  sorption  dehumidifier  program.  The  programs  together,  had 
to  be  interfaced  with  the  use  of  large  common  blocks.  The  grain 
temperature  boundary  of  80F  could  not  be  reasonably  crossed,  due  to  the 
checks  that  would  have  to  be  made  at  each  node  to  determine  which  thin 
layer  equation  to  use  during  transition  periods. 

To  allow  for  variations  in  the  entering  air  properties,  the  ENTRY 
FORTRAN  option  was  used  to  allow  the  calculations  to  begin  inside  the  depth 
loop  after  the  initialization  process  had  taken  place.  Each  hour,  the  air 
exiting  from  the  dehumidifier  had  a  different  temperature  and  humidity.  To 
compensate  for  this  variation,  at  each  hour  the  bottom  nodes  of  the  bed 
were  set  to  the  exit  conditions  of  the  dehumidifier.  (or  the  ambient  air, 
or  collector  exit  conditions,  depending  upon  the  study.)  The  numerical 
analysis  then  continued  in  the  same  manner,  taking  into  account  the  change 
in  the  bottom  node.  These  changes  eventually  were  noticed  at  nodes  higher 
up  in  the  bed  as  the  process  continued.  The  reader  will  notice  when 
consulting  Appendix  I,  that  the  relative  humidity  was  not  reset  to  allow 
for  varying  inlet  conditions.  This  was  because  to  do  so,  would  require 
excessive  calculation  and  this  quantity  was  not  of  primary  concern  to  the 
wri  ter. 

Also  considered  in  the  modification  of  the  fixed  bed  program  were  the 
constants  involved  in  the  calculations,  and  how  they  would  be  affected  by 
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Table  4.B 


THE   EFFECTS  OF  CHANGES   IN  THE  ENTERING  ABSOLUTE  HUMIDITY 
ON  THE  FIXBED  FORTRAN  PROGRAM 


.0001 

A8S 

HUMIDITY 

Name 

.0020 
64.49754 

.0060 
64.08662 

.0100 

CON1 

64.69158 

63.68097 

C0N2 

120.29984 

119.93344 

119.16934 

118.41502 

C0N3 

60.51686 

60.42650 

60.23761 

60.05043 

C0N4 

10.37428 

10.37428 

10.37428 

10.37428 

CON5 

38.71001 

37.71001 

37.71001 

37.71001 

C0N6 

121.03372 

120.85300 

120.47522 

120.10086 

C0N7 

1941.19702 

1942.12183 

1944.04980 

1945.95410 

SC0N1 

86.54869 

86.24405 

85.60966 

84.98454 

SC0N2 

0.44726 

0.44884 

0.45217 

0.45549 

SC0N3 

39.18451 

39.10738 

38.94649 

38.78748 

DELT 

0.02158 

0.02157 

0.02155 

0.02153 

The  change    in  value  o-f   all    constants   is   less   than   Z/.  when    the  entering 
absolute   humidity   wanes   from   0.0001    to  0.0100. 
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changes  in  the  entering  air  properties.  There  are  many  terras  which  do  not 
change  in  a  particular  drying  problem.  The  reader  should  consult  Grain 
Dryer  Simulation  <5)  -for  further  details  concerning  these  constants.  They 
are,  brie-fly,  intermediate  values  that  are  calculated  using  constant  values 
such  as  air  flow  rate,  specific  heats,  and  densities.  It  was  determined 
that  all  constants  changed  in  value  less  than  2%,  when  the  entering 
absolute  humidity  varied  from  0.0001  to  0.0100.  Only  one  constant,  cons7, 
was  found  to  alter  noticeably  with  changes  in  air  temperature.  This  was 
because  cons7  contains  a  &T  term,  and  consequently  varied  considerably  with 
changing  air  temperatures.  Consequently,  all  the  constants  except  cons?, 
were  calculated  only  once.  They  were  calculated  at  the  air  conditions  of 
the  first  hour  of  the  study,  and  remained  constant  throughout.  Cons7, 
however,  was  placed  inside  the  depth  loop,  and  calculated  each  hour  to 
compensate  for  changes  in  entering  air  temperatures,   (see  Table  4.B). 

The  size  of  the  increments  for  depth  and  time  are  extremely  critical 
in  the  fixed  bed  program.  If  too  large,  the  equations  will  diverge  or 
oscillate  from  the  true  solution.  If  too  small,  the  solution  requires 
excessive,  and  expensive  computer  time.  Figure  4.3  contains  approximate 
values  of  air  flow  rate  verses  number  of  layers  per  foot  necessary  for 
stability.  If  the  program  fails  with  good  data,  there  are  two  options 
available  to  the  user;  raise  the  value  of  the  number  of  layers  per  foot,  or 
lower  the  safety  factor  in  the  time  increment  equation. 

The  fixed  bed  program  will  terminate  and  return  control  to  the  main 
program  when  one  of  two  things  happen;  1.)  The  specified  drying  time  has 
been  met.  2.)  The  average  moisture  content  in  the  bed  falls  below  the 
specified  value.  A  sample  output  of  the  fixed  bed  program  as  well  as  the 
modified  listing  can  be  found  in  Appendices  I  and  G,  respectively. 
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DESCRIPTION  OF  SYSTEMS  STUDIED 

Six  different  systems  were  considered  in  the  application  of  solar 
powered  sorption  dehumi  di  f  i  ed  air  to  dry  shelled  corn.  A  brief  outline  of 
each  system  -follows  in  the  next  -few  pages.  Included  in  the  descriptions 
are  schematics  of  the  processes  detailed. 

SYSTEM  1  AND  IA 

The  systems  1  and  IA  are  made  up  o-f  a  fixed  bed  dryer,  that  will  dry 
the  commoditiy,  shelled  corn.  A  fan  is  used  to  blow  air  into  the  plenum 
chamber  located  at  the  bottom  of  the  dryer.  The  air  passes  through  a 
perforated  floor,  then  upward  through  the  grain,  removing  moisture,  and 
exiting  in  a  saturated  state  at  the  top  of  the  dryer.  The  only  difference 
between  the  two  systems  is  the  use  of  an  auxiliary  heater  in  system  IA  to 
preheat  the  inlet  ambient  air,  to  a  set  value  of  100  F.  Modeling  these  two 
systems  will  accomplish  two  things.  First,  the  rate  at  which  the  location 
ambient  air  will  dry  the  grain  can  be  determined.  Secondly,  a  comparison 
can  be  made  between  the  two  systems  on  how  much  faster  the  auxiliary  heater 
will  dry  the  same  bin  of  grain,  and  how  much  the  user  would  have  to  pay  in 
the  way  of  auxiliary  energy  for  that  faster  drying  time.  The  auxiliary 
heat  is  supplied  in  this  study  by  liquefied  petroleum  gas.  The  cost  of 
this  fossil  fuel  was  set  at  *10.00/10*«<S  BTU.  Figures  5.1  and  5.2  are 
schemtics  of  systems  1  and  IA,  respectively. 

The  numerical  modeling  of  the  systems  1  and  IA  required  the  following 
information  and  subroutines  to  be  calculated: 

1.  The  TMY-SOLMET  weather  data  for  the  particular  location 
studied.  (Omaha,  Nebraska) 

2.  The  fixed  bed  dryer  simulation.   This  includes  the  calculation 
of  air  temperature,  product  temperature,  moisture  content, 
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and  absolute  humidity  to  be  calculation  for  each  position 
in  the  bed. 
3.  in  the  case  of  system  1A,  the  calculation  of  the  auxiliary 
energy  required  to  maintain  the  inlet  dryer  air  at  100  F 
dry  bulb  temperature. 

SYSTEM  2  AND  2A 
Systems  2  and  2A  are  represented  in  figures  5.3  and  5.4,  respectively. 
Ambient  air  enters  a  flat  plate  collector  and  is  warmed  to  a  temperature 
dictated  by  the  ambient  conditions  and  collector  variables.  The  warmed  air 
is  then  used  to  dry  the  shelled  corn.  Once  again,  the  difference  between 
the  two  systems  is  that  system  2A  uses  an  auxiliary  heater  to  heat  the  air 
to  a  set  minimum  inlet  temperature  of  100F.  There  will  be  times  when  the 
collector  will  warm  the  air  to  a  temperature  which  is  greater  than  100F. 
If  this  situation  occurs,  the  extra  energy  was  not  allowed  to  be  used.  The 
inlet  air  temperature  remained  constant  at  100F.  This  allowed  the  amount 
of  auxiliary  energy  needed  to  produce  the  same  drying  air  as  the  auxiliary 
heated  ambient  air,  to  be  quantified.  However,  this  is  definately  a 
variable  that  needs  to  be  studied,  as  suggested  in  in  the  section 
RECOMMENDATIONS  FOR  FURTHER  STUDY. 

The  numerical  modeling  of  these  two  systems  required  the  use  of   the 
following  subroutines: 

1.  The  TMY-SOLMET  weather  data  for  the  particular  location 
studied.  (Omaha,  NE) 

2.  The  fixed  bed  dryer  simulation. 

3.  The  collector  performance  subroutine,  COLPER. 

4.  The  calculation  of  the  auxiliary  energy  used  by  system  2A 
to  maintain  the  inlet  dryer  air  at  100F. 

5.  An  economic  analysis  to  determine  the  life  cycle  costs  of 
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a  -flat  plate  collector. 

SYSTEM  3  AND   3A 

The  final  set  of  systems  studied  employed  the  solar  powered  sorption 
dehumidif ier .  Figures  5.5  and  5.6  represent  the  systems  3  and  3A, 
respectively.  The  difference  between  the  two  systems  lies  only  in  the  use 
of  an  auxiliary  heater  in  system  3A  to  insure  the  dryer  inlet  air  is  at 
least  100F  when  entering  the  dryer.  Auxiliary  heat  was  also  used  in  this 
set  of  systems  to  maintain  a  regeneration  temperature  of  176F.  In  system  3 
the  value  calculated  for  the  use  of  auxiliary  heat  is  that  amount  used  to 
regenerate  the  silica  gel  only.  In  the  case  of  system  3A,  the  auxiliary 
heat  costs  were  broken  down  into  two  parts;  the  regeneration  auxiliary 
heat,  and  the  dryer  inlet  auxiliary  heat.  Once  again,  the  inlet  air  in 
system  3A  was  fixed  at  100F,  and  was  not  allowed  to  go  above  that  value. 
This  reduced  the  number  of  variables  to  be  considered  in  the  study. 

The  modeling  of  systems  3  and  3A  used  the  following  subroutines: 

1.  The  TMY-SOLMET  weather  data  for  the  particular  location 
studied. 

2.  The  fixed  bed  dryer  simulation. 

3.  The  collector  performance  subroutine. 

4.  An  economic  analysis  of  the  life  cycle  costs  of  a  solar 
powered  sorption  dehumidifier  compared  to  a  conventional 
dehumidf ier. 

5.  The  simulation  of  a  silica  gel  rotating  bed. 

6.  The  calculation  of  the  auxiliary  energy  costs  needed 
for  the  regeneration  of  the  desiccant,  and  to  maintain 
a  dryer  inlet  temperature  of  100F. 
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RESULTS  AND  DISCUSSION 

To  summarize,  the  six  grain  drying  systems  studied  used  the  -following 
types  of  drying  air: 

1 .  ambient  air 

1A.  auxiliary  heated  ambient  air 

2.  solar  heated  air 

2A.  solar  and  auxiliary  heated  air 

3.  sorption  dehumidified  air 

3A.  sorption  dehumidified  and  auxiliary  heated  air 
Difficulties  arose  in  the  modeling  of  system  2,  the  solar  heated 
drying  air.  This  was  due  to  the  fact  that  the  temperature  of  the  grain 
crossed  the  temperature  boundary  of  80F  twice  within  the  study.  As  was 
discussed  in  the  section  REVIEW  OF  GRAIN  DRYING  MODELS,  different  thiniayer 
equations  would  have  to  be  used  to  correctly  model  the  dryer.  This  would 
lead  to  complicated  checks  at  each  time  increment  and  at  each  node  within 
the  bed  to  specify  which  equation  to  use.  Therefore,  system  2,  the  solar 
heated  air  system,  could  not  reasonably  be  modeled  using  the  present  dryer 
s i mu I  a  t  i  on . 

The  other  five  systems  were  modeled  and  favorable  weather  conditions 
were  sought  in  October  for  Omaha,  Nebraska.  The  TMY-SOLMET  tape  was  read 
for  the  month  of  October  for  Omaha.  The  data  was  reviewed  to  find  a  time 
span  of  about  one  week  which  had  an  average  clearness  index.  The  clearness 
index  is  defined  as  the  ratio  of  the  average  radiation  on  a  horizontal 
surface  to  the  extraterrestrial  radiation.  The  days  selected  for  this 
study  were  October  the  20th  through  the  23rd.  The  parameters  of  the 
systems  were  set  at  the  following  values: 
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1.  regeneration  and  process  air  -flow  rates    120  CFM/FT**2 

2.  density  of  the  air  stream  0.075  lb/FT**3 

3.  rotating  speed  o-f  desiccant  wheel  600  sec/rev 

4.  bed  thickness  0.112  FT 

5.  -frontal  area  0.6092  FT**2 

6.  void  -fraction  0.750 

7.  pressure  drop,  process  &  regeneration  side   4.216  in.  water 

8.  peripheral  leakage  rates  0.0 

9.  collector  area  60.3532  FT**2 

10.  slope  of  collector  10.49  degree 

11.  latitude  41.37  degree 

12.  longitude  96.00  degeee 

13.  duct  temp  drop  to  collector  0.0  F 

14.  duct  temp  drop  to  dehumidif ier  1.8  F 

15.  minimum  temperature  for  regeneration  176.0   F 

16.  heat  exchanger  efficiency  0.0 

17.  income  producing  property 

18.  initial  product  temperature  80.0  F 

19.  bed  depth  5.0  FT 

20.  cross-sectional  area  of  bed  1.0  FT**2 

21.  initial  moisture  content  of  grain  (d.b.)    33.0  '/. 

22.  number  of  layers  per  foot  15.0 

23.  total  time  of  study  96.0   hours 

24.  final  moisture  content  (d.b.)  8.0  '/. 

25.  grain  shelled  corn 

26.  auxiliary  energy  LPG 

27.  cost  of  auxiliary  energy  *1Q. 00/10**6  BTU 

28.  market  discount  rate  12. OX 
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29.  terra  of  mortgage  5  years 

30.  dowpayment  10.  OX 

31.  mortgage  interest  rate  12. OX 

32.  term  of  depreciation  5  years 

33.  property  tax  $0.0 

34.  assessed  valuation  $0.0 

35.  resale  value  *0.0 

36.  miscellaneous  costs  2. OX 

37.  -fixed  equipment  costs  -  solar  41000.00 

38.  -fixed  equipment  costs  -  conventional     $2000.00 

39.  cost  o-f  auxiliary  heater  $250.00 

40.  cost  per  collector  area  $300.00 

41.  -federal  tax  credit  55. OX 

42.  -fuel  inflation  rate  13. OX 

43.  general  inflation  rate  6. OX 

44.  federal  tax  rate  30. OX 

45.  state  tax  rate  5. OX 

46.  term  o-f  economic  analysis  15  years 

47.  accelerated  cost  recovery  depreciation 

The  drying  rates  o-f  the  different  systems  can  be  found  in  Graph  #2. 
The  ambient  air  conditions  for  October,  as  recorded  on  the  TMY-SOLMET 
weather  tape,  will  not  dry  the  corn.  This  will  not  be  true  for  every 
October  in  Omaha,  Nebraska.  However,  for  the  conditions  presented  on  the 
tape,  the  ambient  air  will  not  dry  the  shelled  corn.  In  fact,  the  moisture 
content  of  the  grain  actually  increased  slightly  because  the  EMC  of  the  air 
was  greater  than  the  initial  moisture  content  of  the  corn.  The  other  four 
systems  dried  at  the  rates  shown  in  Graph  tt2.  The  reader  will  notice  that 
the  dehumidified/auxiliary  heated  air  dried  the  fastest.   The  auxiliary 
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heated  ambient  air  UA)  and  the  solar/auxiliary  heated  air  (2A)  dried  at 
the  slowest  rates.  These  two  systems  dried  at  the  same  rates  because  the 
air  entering  the  dryer  was  set  at  100F.  From  this  data  the  savings  in 
auxiliary  energy  due  to  the  use  of  the  solar  collector  can  be  calculated. 
A  conclusion  can  also  be  made  that  the  ambient  air  heated  to  100F  can  not 
dry  as  quickly  as  ambient  air  that  is  dehumidified.  There-fore,  there  is  a 
definite  advantage  to  using  dehumidi-f  i  ed  air  over  heated  air  when 
considering  drying  rates.  The  time  to  dry  the  four  bushels  o-f  corn  <1 
bushel  =  1.25  cu.  ft.)  to  8'/.  moisture  content  (d.b.)  was  much  shorter  for 
the  dehumidified  air.  Both  of  the  systems  that  used  dehumidified  air,  (3> 
and  <3A),  dried  at  about  the  same  rate;  The  dehumidified  and  auxiliary 
heated  air  (3A),  finishing  about  one  hour  before  the  dehumidified  only  air 
(3).  It  should  be  noted  that  an  initial  moisture  content  of  33X  d.b.  and  a 
drying  air  flow  rate  of  120  CFM  is  an  extreme  drying  condition.  Certainly, 
when  applying  this  model  to  predict  the  performance  of  real  dryers,  more 
realistic  parameters  should  be  used. 

A  summary  of  the  auxiliary  energy  usage  can  be  found  in  Figures  6.1, 
6.2,  6.3,  and  6.4.  The  figures  represent  the  daily  costs  for  October  20, 
21,  22,  and  23.  Values  for  the  cost  per  bushel  are  also  presented.  On  the 
first  day,  October  20,  the  auxiliary  energy  costs  are  represented  in  Figure 
6.1.  The  dehumidified/auxiliary  heated  drying  air  required  the  most 
auxiliary  energy  (*3.70  or  *0.?3/bu.).  The  reader  will  notice  that  most  of 
that  energy  was  used  to  regenerate  the  silica  gel.  Only  a  small  portion  of 
the  energy  was  used  to  heat  the  dryer  inlet  air.  Generally,  it  can  be  said 
that  regenerating  the  desiccant  requires  a  lot  of  auxiliary  energy  because 
the  solar  collector  does  not  reach  a  temperature  high  enough  for  sufficient 
regeneration.  This  is  particularly  true  for  the  time  of  year  considered, 
October,  which   is  the  usual  standard  corn  harvesting  time.   Because   the 


FIGURE  6.1 
AUXILIARY  ENERGY  COSTS  FOR  OCTOBER  20 
OMAHA,  NEBRASKA 
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FIGURE  6.2 
AUXILIARY  ENERGY  COSTS  FOR  OCTOBER  21 
OMAHA,  NEBRASKA 
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FIGURE  6.3 
AUXILIARY  ENERGY  COSTS  FOR  OCTOBER  22 
OMAHA,  NEBRASKA 
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FIGURE  6.4 

AUXILIARY  ENERGY  COSTS  FOR  OCTOBER  23 
OMAHA,  NEBRASKA 
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FIGURE   6.5 

TOTAL   COST   OF  AUXILIARY   ENERGY   FOR  THE   STUDY 
OCTOBER   20,21,22,23 
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regeneration  process  requires  a  lot  of  energy,  the  dehumidified  drying  air 
system  required  the  second  highest  amount  of  LPG  <*3.31  or  <*0.83/bu.). 
The  solar/auxiliary  heated  drying  air  required  less  auxiliary  energy  <*1.16 
or  *0.29/bu.)  than  the  heated  ambient  air  <*1.51  or  *0.38/bu.>.  This  is 
expected  because  the  solar  collector  contributes  a  portion  of  the  energy 
needed.  The  ambient  drying  air  did  not  require  an  outside  energy  source, 
but  as  was  -found  in  Graph  tt2,  has  negligible  drying  capacity. 

In  Figure  6.2  the  reader  will  notice  that  the  costs  for  the 
dehumidified  air  systems,  <3>  and  <3A>,  dropped  off  dramatically.  This  is 
because  the  four  bushels  of  corn  being  dried  in  these  two  systems  peached 
their  desired  ending  moisture  content  of  8Z  (d.b.)  about  eight  hours  into 
the  second  day.  As  a  result,  auxiliary  energy  was  not  required  by  these 
two  systems,  <3)  and  <3A),  for  the  remaining  16  hours  of  that  day.  The 
solar/auxiliary  heated  air  and  the  heated  ambient  air,  drying  at  the  same 
rate,  continued  to  use  auxiliary  energy.  The  values  vary  slightly  from  the 
previous  day  due  to  different  weather  conditions  as  supplied  on  the  TMY- 
SOLMET  tape. 

By  the  third  day,  Figure  6.3,  the  dehumidif i ed/auxi 1 iary  heated  system 
(3A)  and  the  dehumidified  system  (3)  did  not  require  any  further  auxiliary 
energy.  Thus,  the  values  for  these  two  systems  are  shown  as  0.0.  Again, 
the  solar/auxiliary  heated  air  and  the  heated  ambient  air  continue  to  dry, 
using  about  the  same  amount  of  energy  as  was  required  in  the  previous  days. 
The  slight  variations  are  due  only  to  changing  weather  conditions.  Figure 
6.4  presents  the  auxiliary  energy  usage  for  the  final  day,  October  23.  The 
same  results  are  obtained  here  as  in  Figure  6.3. 

The  total  costs  of  energy  for  the  study  can  be  found  in  Figure  6.5. 
The  cost  of  drying  per  bushel  is  also  shown.  Because  the  dehumidified  air 
systems,  <3>  and  <3A),  finished  drying  shortly  after  the  second  day  began, 
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FIGURE  6.6 
LIFE  CYCLE  COSTS 
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the  total  amount  of  auxiliary  heat  used  by  the  two  systems  was  actually 
less  overall  than  the  conventionally  heated  ambient  drying  air.  It  should 
be  noted  that  the  pressure  drops  in  the  -fixed  bed  dryer  were  not  considered 
in  this  study.  The  addition  of  this  -factor  would  produce  even  more 
-favorable  results  for  the  sorption  dehumidi-f  ier .  Looking  at  Figure  6.5, 
the  ambient  drying  did  not  use  any  auxiliary  energy,  but  did  not  have  the 
capacity  to  dry  the  grain.  It  is  interesting  to  note  that  the 
solar/auxiliary  heated  drying  air  used  the  least  amount  o-f  auxiliary  energy 
o-f  the  -four  systems  that  used  it.  There  is  a  different  of  about  10  cents 
per  bushel  to  dry  the  corn  using  the  dehumidified/heated  air  versus  the 
dehumidified  only  air.  This  differece  represents  the  cost  to  heat  the 
drying  air.  The  rest,  which  is  about  9\V.  of  the  cost,  is  due  to  the 
auxiliary  heat  needed  to  regenerate  the  silica  gel.  If  this  number  could 
somehow  be  lowered,  this  method  of  drying  grain  would  become  more 
attractive. 

When  comparing  drying  rates  and  their  costs,  it  is  apparent  that  the 
fastest  drying  rate  is  obtained  with  the  dehumidified/auxiliary  heated  air. 
In  Graph  *2,  the  reader  will  notice  that  the  dehumidified/auxiliary  heated 
and  the  dehumidifed  drying  air  reach  the  desired  ending  moisture  content  of 
&/.  d.b.  at  about  the  same  time.  There  was  infact,  only  a  difference  of  one 
hour  between  the  two  drying  times.  The  dehumidified  air  dries  at  a  faster 
rate  and  is  also  less  expensive  than  the  conventionally  heated  ambient  air. 
However,  the  equipment  needed  to  use  a  solar  powered  sorption  dehumidifier 
is  expensive.  To  fairly  determine  the  feasibility  of  the  system  will 
require  an  economic  analysis. 

The  last  phase  of  this  study  compared  the  life  cycle  costs  of  the 
systems.  This  value  is  the  sum  of  all  the  costs  associated  with  the  energy 
delivery  system  over  its  lifetime  or  over  a  selected  period  of  analysis,  in 
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today's  dollars,  and  takes  into  account  the  time  value  of  money.  The  basic 
idea  o-f  life  cycle  costs  is  that  anticipated  future  costs  are  brought  back 
to  present  cost  (discounted)  by  calculating  how  much  would  have  to  be 
invested  at  a  market  discount  rate  to  have  the  -funds  available  when  they 
will  be  needed.  The  market  discount  rate  is  the  pate  of  return  on  the  best 
alternative  investment.  The  dryer,  fan,  and  duct  work  were  considered 
common  to  all  the  systems.  Therefore,  the  numbers  presented  in  Figure  i.6 
represent  the  difference  between  the  base  unit  and  the  addition  of  the 
solar  related  equipment. 

The  solar  powered  sorption  dehumidifier  with  the  auxiliary  heater  had 
the  largest  solar  life  cycle  cost  (SLCO  of  $2344.  The  sorption 
dehumidifier  without  the  auxiliary  heater  had  a  SLCC  of  $2150  for  the  15 
year  economic  analysis.  The  solar  collector  and  auxiliary  heater  had  a 
SLCC  of  $558.  The  heated  ambient  air,  which  required  the  purchase  of  an 
auxiliary  heater,  had  a  conventional  life  cycle  cost  (CLCO  of  $200.00. 
The  ambient  air  system  has  a  *0.0  CLCC  since  the  equipment  used  is  that 
equipment  common  to  all  the  systems. 

Are  there  any  circumstances  under  which  the  sorption  dehumidifier 
might  be  economically  feasible?  The  restriction  of  running  the  solar 
energy  powered  dehumidifier  for  the  four  days  of  this  study  makes  it  very 
difficult  for  the  investment  to  be  attractive.  Even  the  use  of  the 
sorption  dehumidifier  for  the  two  to  four  weeks  of  corn  harvest  that  occur 
each  year  would  doubtfully  produce  economically  feasible  conditions.  The 
need  for  other  on-farm  applications  is  great.  Some  examples  of  other  uses 
for  the  solar  energy  powered  sorption  dehumidifier  would  include: 

1.)  The  continuous  regeneration  of  silica  gel  beds  that 
could  be  stored  until  needed. 

2.)  If  the  dehumidifier  could  be  considered  by  components, 
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then  certainly  the  solar  collector  alone,  could 
■find  use  in  supplying  warm  air  for  animal  shelters. 
The  solar  energy  powered  sorption  dehumidifier  computer  model  is  now 
in  a  -form  that  can  be  used  to  investigate  many  parameter  variations.   For 
instance,   how  most  energy  would  have  to  be  supplied  to  ambient  air  in  the 
form  of  heat,  to  get  equivalent  drying  rates  between  heated  and  dehumidifed 
air?  The  number  o-f  studies  that  could  be  performed  are  numerous. 

The  GRAIN  DRYER  SIMULATION  coupled  with  the  dehumidifier  model,  the 
collector  model,  the  economic  analysis,  and  the  weather  data  information 
require  about  30  minutes  (depending  upon  the  specifics  of  the  study)  of 
computer  time.  Of  the  average  30  minutes  required,  approximately  50%  of 
that  time  was  spent  in  the  fixed  bed  model  dryer.  This  study  could  not 
include  a  more  indepth  analysis  due  to  the  lack  of  a  funding  source  for  the 
expensive  computer  time  required  on  the  University's  large  computer. 
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SUMMARY  AND  CONCLUSIONS 

The  main  objective  of  this  study  was  to  determine  a  method  of 
simulating  grain  drying  with  a  solar  powered  sorption  dehumi di tier.  The 
type  of  grain  dryer  selected  was  -fixed  bed,  and  the  model  used  was  the 
GRAIN  DRYER  SIMULATION  developed  at  Michigan  State  University  by  Bakker- 
Arkema,  et  al .  Limitations  of  this  model  include: 

1.)  The  product  temperature  can  not  cross  the 

temperature  boundary  of  80F  without  resulting 
added  computer  time  and  expense,  due  to  the 
multitude  of  checks  that  would  have  to  be  made 
to  be  certain  the  correct  equation  was  being  used 
at  each  node  of  the  bed. 
2.)  The  dryer  simulation  requires  a  lot  of  time  on  the 
the  National  Advanced  Systems/6620  at  KSU. 
For  the  four  day  study,  the  time  required  on  the 
ccmputer  ranged  from  15  to  40  minutes,  depending 
on  the  system  modeled.  Some  of  this  time  was 
spent  in  modeling  the  dehumidifer  bed,  the  collector 
and  reading  the  weather  data.   But,  over  5QY.   of  the 
total  time  was  required  by  the  dryer  model.   While 
this  model  is  good,  it  is  not  practical  in  its  current 
state.  The  writer  has  made  a  number  of  suggestions  to 
remedy  this  problem  in  RECOMMENDATIONS  FOR  FURTHER 
STUDY . 
3.)  This  grain  dryer  simulation  was  developed  to  predict 
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the  drying  characteristics  of  a  certain  type  of  dryer 
when  the  entering  air  properties,  such  as  temperature 
and  humidity,  remained  constant.   The  changing  of 
the  inlet  air  properites  on  an  hourly  basis,  and  how 
this  effects  the  validity  of  the  model  are  unknown. 
Research  wi 1 1  haue  to  be  performed  to  compare 
experimental  results  with  the  simulated  results. 
Within   the  context  of   this  study,   the  solar  powered  sorption 
dehumidifier  required  an  additional   investment   of   about   52144.00,   over 
convent ional  ,   fossil  fuel  dried  grain.   If  this  system  is  to  be  considered 
feasible,  other  on-farm  applications  will  haue  to  be  found  in   addition   to 
the  drying  of  grain. 
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RECOMMENDATIONS  FOR  FURTHER  STUDY 

1.)  Due  to  the  expense  o-f  running  this  program  on  the  National  Advanced 
Systems/6620  at  K'SU,  one  o-f  two  suggestions  are  made: 

a.  Decide  upon  a  microcomputer  system  and  trans-fer  the 
program  to  that  system  so  that  the  money  -factor  is 
essentially  eliminated.   Access  to  weather  data  will 
still  be  needed  and  a  method  will  have  to  be  determined 
to  obtain  that  data  -  whether  by  the  use  o-f  a  modem, 

or  putting  the  data  on  a  -floppy  disk. 

b.  Consider  a  more  simpli-fied  deep  bed  drying  model  that 
is  empirical  in  nature  and  can  be  applied  to  drying 
shelled  corn,   <However,  this  will  limit  the  application 
to  corn  drying  unless  individual  empirical  models  can  be 
•found  -for  each  drying  use.) 

2.)  Once  the  computer  model  is  in  a  less  expensive  -form,  obtain  the 
equipment  and  instrumentation  needed  to  validate  the  numerical 
mode  1 . 

3.)  Model  Concurrent,  Cross-flow,  and  Counter-flow  dryers,  with  the  sorption 
dehumi  di -f  i  ed  air,  and  validate  the  models  with  the  necessary  equipment 
and  instrumentation. 

4.)  Develop  a  design  computer  model  to  size  the  Solar  Energy  Powered 
Sorption  Dehumi  di  -f  i  cat  i  on  System. 
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APPENDIX  A 
OPTIMUM  ANGLE  FOR  A  FLAT  PLATE  SOLAR  COLLECTOR  FACING  DUE  SOUTH 


FOR  A  FLAT  PLATE  COLLECTOR  SLOPED  TO  THE  SOUTH,  THE  ANGLE  OF  INCIDENCE  IS 
THE  ANGLE  BETWEEN  THE  BEAM  RADIATION  ON  THE  SURFACE  OF  THE  COLLECTOR,  AND 
THE  NORMAL  TO  THAT  SURFACE. 


6  =  ANGLE  OF  INCIDENCE 


FOR  LOCATIONS  IN  THE  NORTHERN  HEMISPHERE: 


COS 


6  =  COS(  ^  -  B  )  COS  C  COS  UT  *   SIN<  (f>    -  8  >  SIN  C 

WHERE: 

Q    =  ANGLE  OF  INCIDENCE 

ft    =  LATITUDE 

5  =  DECLINATION 

p  =  SLOPE 

W  =  HOUR  ANGLE 


DECLINATION 


DECLINATION   IS  THE  ANGULAR  POSITION  OF  THE  SUN  AT  SOLAR  NOON  WITH  RESPECT 
TO  THE  PLANE  OF  THE  EQUATOR,  NORTH  POSIT1UE.    <-23.45  <  C.  <  23.45) 


=23.45  SIN 


3,60 


2fl4_i_fl 
365 


NUMBER   OF  THE   DAY   OF  THE  YEAR 


HOUR  ANGLE 

THE  HOUR  ANGLE  IS  THE  ANGULAR  DISPLACEMENT  OF  THE  SUN  EAST  OR  WEST  OF  THE 
LOCAL  MERIDIAN  DUE  TO  ROTATION  OF  THE  EARTH  ON  ITS  AXIS  AT  15  DEGREES  PER 
HOUR,  MORNING  NEGATIVE,  AFTERNOON  POSITIVE. 
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MAXIMIZE  INCIDENT  RADIATION 

TO  MAXIMIZE  THE  INCIDENT  RADIATION,  MINIMIZE  6    (MAXIMIZE  COS  8   ). 


dCQS-A =  -SINC0-  B  >C03  Cj  COS  W  *      COS<  (j)  -  B  >SIN  C.  =0 


?-p> 

=  TAN 

S- 

cos  ur 

-  0- 

_! 
•  TAN 

_  cos  or 

REARRANGING: 

TAN( 


FOR  OPTIMUM  SLOPE  USE  THE  DECLINATION  ANGLE  THAT  CORRESPONDS  TO  THE  AVERAGE 
FOR  THE  MONTH  THE  STUDY  COVERS.  (SEE  DUFF IE  AND  BECKMAN,  PAGE  12> 


TO  FIND  THE  MEAN  VALUE  FOR  COS  UT   ,  AN   INTEGRABLE  FUNCTION,   USE  THE 
DEFINITION: 

F    =   Jo.  Eixldx 
(b-  ») 

THEREFORE :  T/4 

COS  U   =2   Jo  CDSi£dML_        =   0.900 
TT/4  -  <-7T/4> 
THIS  VALUE  FOR  THE  MEAN  VALUE  OF  COS  UT  ,  IS  THE  AVERAGE  FOR  EVERYDAY. 

THE  VALUE  THAT  SHOULD  BE  ENTERED  IN  THE  DATA  INPUT  LINE  IS:  SL 


ft 


SL  =   t>    "  LATITUDE 
opr 
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APPENDIX  B 
DeBOER  EMC  EQUATION  FOR  SHELLED  CORN 

The   following  empirical    equations  are  -for    the   EMC  values  of   shelled  corn. 

Me  =(Sl*rh**3>/1.02   ♦    CF1/0.17  -  0 .028333*S! )rh  0.0    <   rh    <   0.17 

Me   =SK0.34-rh)**3/1.02   +   S2<rh-0. 17)**3/1 .02   +    (F2/0 . 17-0 .028333S2)*<rh-0 . 17) 
+    (Fl/0.17  -   0.028333Sl)*(0.34-rh)  0.17   <   rh    <    0.34 

Me  =S2<0.51-rh)*»3/l .02   +   F3/0 . 17<rh-0 .34)    +    (F2/0  .17-0.028333S2WO  .51-rh> 

0.34    <  rh    <    0.50 

Me  =   S3<rh-0.49)**3/l .02   +    (F5/0.17   -  0 ,028333*S3]*<rh-0 .49)    +   F4/0.17* 
f.0. 66   -rh)  0.50    <   rh    <    0.66 

Me  =  S3(0.83-rh)**3/l  .02  +  S4(rh-0  .66)**3/l  .02  +tF6/0.17  -0  .028333S41* 
(rh-0.64)  +  [F5/0.17  -  0 .028333S31*(0 .83-rh)    0.66   <  rh  <  0.83 

Me  =  S4(l .00-rh)**3/l .02  +  F7/0 . 17<rh-0 .83)  +[F6/0.17  -  0 .028333S4]*( 1 .00-rh) 

0.83  <  rh  <  1 .00 

where : 

rh  =  relative  humidity,   decimal 

Fl   =   -0.0003922T   +   0.1000 

F2  =   -0.0004353T   +   0.1328 

F3  =   -0.0005359T   +   0.164<S 

F4  =   -0.0005375T   +    0.1624 

F5  =   -0.0007075T   +   0.2075 

F6  =   -0.0007449T   +   0.2532 

F7  =   -0.001071T     +   0.3931 

51  =   13.33(-9F1    t   6F2   -  F3) 

52  =   13.83(4F3   -   9F2   4   6F1 ) 

53  =    13.33(4F4   -   9F5  +    6F6  -  F7) 

54  =   13.83<4F7  -   9F6   +   4F5   -  F4) 
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APPENDIX  C 

THOMPSON  EQWTION  FOR  EMC 

The  Thompson  Equation  -for  determining  EMC,  used  in  the  fixed  bed  program 
for  temperatures  above  235F. 

EMC  =  0.01*SQRT<(-ALOG<1.0-RH))/(0.0000382*(T+50.0J)) 


WHERE 

EMC  =  EQUILIBRIUM  MOISTURE  CONTENT 
RH  =  RELATIVE  HUMIDITY,  DECIMAL 
T   =  TEMPERATURE  OF  THE  AIR 


16 


APPENDIX  D 
THE  SUBBAH  EMPIRICAL   DRYING  EQUATION  FOR   CORN 

TEMPERATURE   RANGE:      3*   -   70F 

MR  =   EXP   [    -k<    t**0.4e4>] 

where 

k     =      exp    (-x*t**y) 

x     =      [6.0142+    (1  .453E-04Xrh>**2]**0.5 

-   «   l«t3.353E-M   +    <3.0E-08>(rh>**2]**0.5 

y     =      0. 1245-2. 197E-03(rh>    +   2.3E-05<rh>e-5.8E-05*0 
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APPENDIX  E 

SHELLED  CORN  DRYING  EQUATIONS  FOR  THE  TEMPERATURE  RANGE  90  -  160F 
TROEGER  AND  HUKILL 


t/60   =   pKMbar  -  Me)**ql  -  pHMo  -  Me)**ql  Mo  >  Mbar  >Mxl 

1/60  =   p2<Mbar  -  Me)**q2  -  p2<Mxl  -  Me)**q2  *    txl       Mxl  >  Mbar  >  Me 
i/60   =  p3<Mbar  -  Me)**q3  -  p3<Mx2  -  Me)**q3  +  tx2       Mx2  >  Mbar  >  Me 

WHERE 

Mxl  =  0.40<Min  -  Me)  +  Me 

Mx2  =  0.12(Min  -  Me)  +  Me 

txl  =  [pi  (Mxl  -  Me)**ql  -  pKMin  -  Me)»*ql]/«0 

tx2  =  Cp2<Mx2  -  Me)**q2  -  p2<Mxl  -  Me)**q2]/40  +  txl 

pi   =  expC-2.45  -  (6.42*Min**l .25)  -  3.15*rh  +  (9.i2«Min* 
rh**.5)  +  0.030*6  -  0.12*Ua) 

p2  =  exp[  2.82  +  7.49<rh  +0 .01 >**0 .67   -  0.01799] 

p3  =  0.12C (Min  -  Me)**Cq2-q3) I*<p2»q2/q3) 

ql  =-3.98  +  2.37Min  -  [0.019/(rh  +  0.015)1  +  0.01*9 

q2  =-exp<0.810  -  3. 1 lrh) 

q3  =  -1 .0 
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APPENDIX  F 

COMPUTER  LISTING  OF  THE  MAIN  PROGRAM  USED  IN  THE  MODELING  OF 

THE  SOLAR  POWERED  SORPTION  DEHUMIDIFICATION  SYSTEM 

The  main  program  is  listed  here  due  to  a  number  of  changes  which  were 
made  during  the  writer's  study.   The  changes  include: 

1.  The  inclusion  of  new  variables  in  the  common  data  blocks 
to  allow  inter-facing  with  the  fi:<ed  bed  dryer  simulation. 

2.  Coding  to  allow  for  the  calculation  and  output  of  the 
auxiliary  costs  used  to  heat  the  dryer  inlet  air,  and  to 
regenerate  the  desiccant.   This  was  performed  on  a  daily 
basis  and  a  total  amount  was  also  calculated  for  the  study. 

3.  Alterations  were  made  to  the  algorithm  needed  to  read  the 
TMY-SOLMET  Weather  Tape.   The  program  now  has  the  capabilities 
to  read  across  month  boundaries  and  to  read  the  first  day  of 
the  year. 

4.  Using  the  amount  of  moisture  removed  during  the  drying  process, 
the  data  from  the  dryer  simulation  was  used  to  calculate  the 
fraction  of  the  load  supplied  by  the  system  configuration. 

5.  Coding  to  enter  the  dryer  simulation,  fixbed,  and  to 

enter  the  depth  loop,  sub-fix,  during  subsequent  calculations. 
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APPENDIX  G 

THE  FIXED  BED  DRYER  SIMULATION 
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APPENDIX  H 


GLOSSARY  FOR  MAIN  VARIABLES  USED  IN  SYSTEMS  1,  1A,  2,  2A,  3,  AND  3A 

AC  =  AREA  OF  COLLECTOR 

AF  =  FRONTAL  AREA  OF  BED 

ALONG  =  COLLECTOR  LOCATION  LONGITUDE 

AMASSP  =  AIR  MASS  FLOU  RATE,  PROCESS  SIDE 

AMASSR  =  AIR  MASS  FLOW  RATE,  REGENERATION  SIDE 

AY  =  INTERNAL  SURFACE  AREA  OF  DESSICANT  PER  UNIT  VOLUME 

AVFR  =  AIR  VOLUME  FLOW  RATE  IN  LITERS  PER  SECOND 

BETA  =  INTERMEDIATE  VALUE  USED  TO  DETERMINE  AMBIENT  ABSOLUTE  HUMIDITY 

BETAF  =  FRACTIONAL  FORM  OF  BETA 

BLEAK  =  BYPASS  LEAKAGE  RATE 

BTUH20  »  ENERGY  INPUT  PER  WATER  REMOVED  (BTU/LB-H20) 

CA  =  UNIT  COST  VARYING  WITH  SIZE  OF  SOLAR  ENERGY  SYSTEM 

CAP  =  HEAT  CAPACITY  OF  DRY  PRODUCT 

CAR  =  HEAT  CAPACITY  OF  DRY  AIR 

CFM  =  AIRFLOW  RATE  AT  INLET  AIR  TEMPERATURE  FT**3/MIN/FT**2 

CEC  =  EQUIPMENT  COST  -  CONVENTIONAL 

CES  =  EQUIPMENT  COST  -  SOLAR 

CF  =  COST  OF  FUEL  IN  DOLLARS 

CFLAG  =  PROGRAMMING  FLAG  TO  DESIGNATE  THE  DIFFERENCE  BETWEEN  INCOME 

PRODUCING  AND  NON-INCOME  PRODUCING  SOLAR  SYSTEMS 

CLK  =  CLOCK  HOUR  UNDER  CONSIDERATION 

CMC  =  CONVENTIONAL  RATIO  OF  MISC.  COSTS  TO  INITIAL  INVESTMENT 

CMCP  =  PERCENTAGE  FORM  OF  CMC 

CMIR  =  ANNUAL  MORTGAGE  INTEREST  RATE  -  CONVENTIONAL 

CMIRP  =  PERCENTAGE  FORM  OF  CMIR 

COST  =  DAILY  COST  OF  AUXILIARY  HEAT  SUPPLIED  AT  DRYER  INLET 

CV  =  HEAT  CAPACITY  OF  WATER  VAPOR  (BTU/LB-F) 

CW  =  HEAT  CAPACITY  OF  LIQUID  WATER  (BTU/LB-F) 

DAY  =  DAY  UNDER  CONSIDERATION 

DAYB  =  BEGINNING  DAY  OF  STUDY 

DAYF  =  CONTROL  TO  AID  IN  OBTAINING  THE  PREVIOUS  DAY  NUMBER  WHEN  CROSSING 

MONTH  BOUNDARIES 

DBT  =  DRY  BULB  TEMPERATURE 

DBTPR  =  DEPTH  BETWEEN  DESIRED  OUTPUT  IN  X-DIRECTION 

DC  =  DOWNPAYMENT  -  CONVENTIONAL 

DCP  =  PERCENTAGE  FORM  OF  DC 

DD  =  DESICCANT  CHAR  DIMENSION 

DEEP  =  DEPTHS  AT  WHICH  OUTPUT  OCCURS   (FT) 

DELT  =  TIME  INCREMENT  (HR) 

DELX  =  DEPTH  INCREMENT  OR  WIDTH  INCREMENT  (FT) 

DEND  »  ENDING  DAY  OF  STUDY 

DEPC  =  DEPRECIATION  TAX  DEDUCTION  -  CONVENTIONAL 

DEPS  =  DEPRECIATION  TAX  DEDUCTION  -  SOLAR 

DEPTH  =  TOTAL  BED  DEPTH 

DIFR  =  DIFFUSE  RADIATION 

DIRR  =  DIRECT  RADIATION 

DIC  =  DEPRECIATION  INFLATION  RATE  -  USED  FOR  DECLINING  BALANCE  -  CONVENTIONAL 

DIS  =  DEPRECIATION  INFLATION  RATE  -  USED  FOR  DECLINING  BALANCE  -  SOLAR 

DMCC  =  DISCOUNTED  VALUE  OF  MISC.  COSTS  -  CONVENTIONAL  < INSURANCE. ETC. ) 

DMCS  =  DISCOUNTED  VALUE  OF  MISC.  COSTS  -  SOLAR 

DOPC  =  FLAG  INDICATING  WHICH  DEPRECIATION  SCHEDULE  IS  UTILI ZED-CONVENTIONAL 
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DOPS  =  FLAG  INDICATING  WHICH  DEPRECIATION  SCHEDULE  IS  UTI LI 2 ED-SOLAR 

DPT  =  DEW  POINT  TEMPERATURE 

DOAUX  =  AUXILIARY  HEAT  SUPPLIED  AT  INLET  OF  DRYER 

DRC  =  MARKET  DISCOUNT  RATE  -  CONVENTIONAL 

DRCP  =  PERCENTAGE  FORM  OF  DRC 

DRS  =  MARKET  DISCOUNT  RATE  -  SOLAR 

DRSP  =  PERCENTAGE  FORM  OF  DRS 

DRVC  =  DISCOUNTED  RESALE  VALUE  -  CONVENTIONAL 

DRVS  =  DISCOUNTED  RESALE  VALUE  -  SOLAR 

DS  =  DOUNPAYMENT  -  SOLAR 

DSP  =  PERCENTAGE  FORM  OF  DS 

DTCD  =  TEMP  DROP  THROUGH  DUCTING,  COLLECTOR  TO  DEHUMIDIFIER 

DTDC  =  TEMP  DROP  THROUGH  DUCTING.  DEHUMIDIFIER  TO  COLLECTOR 

DTQU  =  DAILY  TOTAL  HEAT  GAINED  IN  THE  COLLECTOR 

DTQAUX  =  DAILY  TOTAL  AUX  HEAT  SUPPLIED  FOR  REGENERATION  OF  DESICCANT 

EFF  =  EFFICIENCY 

ENCR  =  ENGINEERING  CORRECTED  RADIATION 

ENERGY  =  ENERGY  INPUT  (CUMMULATIVE)       <BTU/HR-FT**2> 

EXTR  =  EXTRATERRESTRIAL  RADIATION 

F  =  FRACTION  OF  HEAT  SUPPLIED  BY  SOLAR,  WITH  RESPECT  TO  THE  TOTAL  HEAT 

SUPPLIED  (1-F  IS  THE  PERCENT  AUXILIARY  AND  ADDITIONAL  HEAT  SUPPLIED! 

FF  =  PERCENTAGE  FORM  OF  F 

FI  =  FUEL  INFLATION  RATE 

FIP  =  PERCENTAGE  FORM  OF  FI 

FL  =  LENGTH  OF  DESICCANT  WHEEL  IN  FLOW  DIRECTION 

FP  =  COLLECTOR  FACTOR  FRUL 

FR  =  COLLECTOR  FACTOR  FRTA 

FRIC  =  DESICCANT  FRICTION  FACTOR 

FTR  =  FEDERAL  TAX  RATE 

FTRP  =  PERCENTAGE  FORM  OF  FTR 

GA  =  DRY  AIRFLOW  RATE 

GEMT  =  GRAND  TOTAL  EVAPORATED  MOISTURE 

GI  =  GENERAL  INFLATION  RATE 

GIP  ■  PERCENTAGE  FORM  OF  GI 

GTQU  =  GRAND  TOTAL  OF  SOLAR  HEAT  GAIN 

GTQAUX  =  GRAND  TOTAL  OF  AUX.  HEAT  USED  FOR  REGENERATION  OF  DESICCANT 

HC  =  CONVECTIVE  HEAT  TRANSFER  COEFFICIENT  <BTU/HR-FT**2-F> 

HFG  =  LATENT  HEAT  OF  WATER  IN  GRAIN  (BTU/LB) 

HIN  =  INLET  HUMIDITY  RATIO   <LB-H20/LB-DRYAIR> 

HRB  =  BEGINNING  HOUR  OF  STUDY 

HEND  =  ENDING  HOUR  OF  STUDY 

HMU  =  AIR  TO  DESICCANT  MASS  RATIO  IN  DESICCANT  WHEEL 

1  =  DAY  COUNTER 

I DATE  =  NUMBER  DAY  OF  THE  YEAR 

IFLAG  =  FLAG  TO  SIGNAL  AUXILIARY  HEATER  INTO  DEHUMIDIFIER 

IHR  =  SOLAR  HOUR 

ILEAP  =  PROGRAM  AID  TO  DETERMINE  THE  NUMBER  DAY  OF  THE  YEAR  TAKING 

INTO  CONSIDERATION  LEAP  YEAR 

INDPR  =  NUMBER  OF  NODES  BETWEEN  PRINTS 

IZ  =  CONTROL  TO  AID  IN  SETTING  TAPE  TO  THE  2ERO  HOUR  OF  THE  DAY  DATA 

IS  TO  BE  READ 

J  =  HOUR  COUNTER 

KT  =  HOURLY  CLEARNESS  INDEX 

LADD  =  ADDITIONAL  LOAD  VARIED  TO  DETERMINE  OPTIMUM  LIFE  CYCLE  COST 

LCCC  =  LIFE  CYCLE  COSTS  -  CONVENTIONAL 

LCCS  =  LIFE  CYCLE  COSTS  -  SOLAR 


LP 

LPP 

MEND 

MHR 

MINS 

MO 

MOB 

MOTEST 

MPIC 

MP1S 

ND 

NDAY 

NDEPC 

NDEPS 

NDEPPC 

NDEPPS 

NEA 

NH 

NLC 

NLPF 

NLS 

NMINC 

NMINS 

NMINPC 

NMINPS  = 


OBSR 
PATM 
PDAY 

PERB 

PERD 

PHI 

PLEAKP 

PLEAKR 

PONEC 

PONES 

PMO 

PSEA 

PSTN 

PSTNC 

PTCC 

PTCS 

PTUOC 

PTUOS 

PW_C 

PWF 

PWFP 

PW_S 


D.  -  CONVENTIONAL 
D.  -  SOLAR 
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=  LOAD  DRIED  BY  DEHUMIDIFIER 

=  LOAD  DRIED  BY  DEHUMIDIFIER  AND  SUPPLEMENTAL  SOURCE 

=  ENDING  MONTH  OF  STUDY 

=  TAPE  HOUR  READ 

=  MINUTES  OF  SUNSHINE 

=  MONTH 

=  BEGINNING  MONTH  OF  STUDY 

=  PROGRAM  AID  TO  DETERMINE  THE  NUMBER  DAY  OF  THE  YEAR 

=  DISCOUNTED  MORTGAGE  PRINCIPAL  AND  INTEREST  -  CONVENTIONAL 

=  DISCOUNTED  MORTGAGE  PRINCIPAL  AND  INTEREST  -  SOLAR 

=  NUMBER  OF  DAYS  INCLUDED  IN  STUDY 

=  NUMBER  OF  DAY  IN  A  GIVEN  MONTH 

=  TERM  OF  DEPRECIATION  -  CONVENTIONAL 

=  TERM  OF  DEPRECIATION  -  SOLAR 

=  TERM  OF  DEPRECIATION,  USED  FOR  D.  B.  AND  S.  0. 

=  TERM  OF  DEPRECIATION,  USED  FOR  D.  B.  AND  S.  0. 

=  TERM  OF  THE  ECONOMIC  ANALYSIS 

=  NUMBER  OF  HOURS  IN  THE  STUDY 

=  TERM  OF  THE  LOAD  -  CONVENTIONAL 

=  NUMBER  OF  LAYERS  PER  FOOT 

=  TERM  OF  THE  LOAN  -  SOLAR 

=  YEARS  THE  MORTGAGE  PAYMENTS  CONTRIBUTE  TO  THE  ANALYSIS  -  CONVENTIONAL 

EQUAL  TO  THE  MINIMUM  OF  NLC  AND  NEA 
=  YEARS  THE  MORTGAGE  PAYMENTS  CONTRIBUTE  TO  THE  ANALYSIS  -  SOLAR 

EQUAL  TO  THE  MINIMUM  OF  NLS  AND  NEA 
=  YEARS  THE  MORTGAGE  PAYMENTS  CONTRIBUTE  TO  THE  ANALYSIS  FOR  THE 

EQUIPMENT  -  CONVENTIONAL 

YEARS  THE  MORTGAGE  PAYMENTS  CONTRIBUTE  TO  THE  ANALYSIS  FOR  THE 

EQUIPMENT  -  SOLAR 

EQUAL  TO  THE  MINIMUM  OF  NEA  AND  NDEPS 

OBSERVED  RADIATION 

ATMOSPHERIC  PRESSURE 

CALENDAR  NUMBER  OF  THE  DAY  PRECEDING  THE  REQUESTED  BEGINNING  DAY 

OF  STUDY 

PERCENTAGE  OF  BEAM  RADIATION  BASED  ON  EXTRATERRESTRIAL  RADIATION 

PERCENTAGE  OF  DIFFUSE  RADIATION  BASED  ON  EXTRATERRESTRIAL  RADIATION 

COLLECTOR  LOCATION  LATITUDE 

PERIPHERY  LEAKAGE  RATE,  PROCESS  SIDE 

PERIPHERY  LEAKAGE  RATE,  REGENERATION  SIDE 

RATIO  OF  LIFE  CYCLE  FUEL  COSTS  TO  FIRST  YEAR  FUEL  COSTS  -  CONVENTIONAL 

RATIO  OF  LIFE  CYCLE  FUEL  COSTS  TO  FIRST  YEAR  FUEL  COSTS  -  SOLAR 

NUMBER  OF  THE  MONTH  PRECEDING  THE  REQUESTED  BEGINNING  MONTH 

ATMOSPHERIC  PRESSURE  AT  SEA  LEVEL 

ATMOSPHERIC  PRESSURE  AT  THE  STATION 

ATMOSPHERIC  PRESSURE  AT  THE  STATION  CORRECTED 

DISCOUNTED  PROPERTY  TAX  COSTS  -  CONVENTIONAL 

DISCOUNTED  PROPERTY  TAX  COSTS  -  SOLAR 

LIFE  CYCLE  COSTS  OF  ADDITIONAL  CAPITAL  INVESTMENT  TO  INITIAL 

INVESTMENT  -  CONVENTIONAL 

LIFE  CYCLE  COSTS  OF  ADDITIONAL  CAPITAL  INVESTMENT  TO  INITIAL 

INVESTMENT  -  SOLAR 

PRESENT  WORTH  FACTOR  FOR  CALCULATING  A  TERM  IN  PONE  OR  PTWO  - 

CONVENTIONAL 

PRESENT  WORTH  FACTOR,  SEE  PG  386  DUFF  IE  AND  BECKMAN 

SECOND  FORM  OF  PWF,  SEE  ABOVE 

PRESENT  WORTH  FACTOR  FOR  CALCULATING  A  TERM  IN  PONE  OR  PTWO  - 


lid 


SOLAR 

PUW  =  PARTIAL  PRESSURE  OF  WATER  UAPOR 

QDEHU  =  RATE  OF  HEAT  TRANSFERRED  IN  DEHUMIDIFIER,  WATTS 

QDUCT  =  RATE  OF  HEAT  LOST  FROM  DUCTING,  WATTS 

QHEXP  =  RATE  OF  HEAT  TRANSFERRED  IN  HEAT  EXCHANGER,  WATTS 

QU  =  RATE  OF  HEAT  GAINED  IN  COLLECTOR.  WATTS 

QAUX  =  RATE  OF  HEAT  SUPPLIED  BY  AUXILIARY  HEATER  INTO  DEHUMIDIFIER,  WATTS 

QSURP  =  SURPLUS  SOLAR  HEAT  AVAILABLE,  WATTS 

REFLEC  =  GROUND  RELECTANCE  RATIO 

RH  =  RELATIVE  HUMIDITY 

RHC  =  SATURATION  RELATIVE  HUMIDITY  =  0.9999999999 

RHIN  =  INLET  RELATIVE  HUMIDITY   (DECIMAL) 

RHOF  =  AVERAGE  DENSITY  OF  AIR  FLOW  STREAM 

RHOP  =  DRY  BULK  DENSITY  OF  GRAIN   (LB/FT**3> 

RLEAK  =  RECIRCULATION  LEAK  RATE 

RPS  =  REV  PER  SEC 

RVC  =  RATIO  OF  RESALE  VALUE  AT  END  OF  PERIOD  OF  ANALYSIS  TO  INITIAL 

INVESTMENT  -  CONVENTIONAL 

RVS  =  RATIO  OF  RESALE  VALUE  AT  END  OF  PERIOD  OF  ANALYSIS  TO  INITIAL 

INVESTMENT  -  SOLAR 

SKY  =  SKY  CONDITIONS,  TYPE  AND  EXTENT  OF  CLOUD  COVER 

SL  =  SLOPE  OF  COLLECTOR  (OPTIMUM  ANGLE-LATITUDE) 

SMC  =  SOLAR  RATIO  OF  MISC.  COSTS  TO  INITIAL  INVESTMENT 

SMCP  =  PERCENTAGE  FORM  OF  SMC 

SMIR  =  ANNUAL  MORTGAGE  INTEREST  RATE  -  SOLAR 

SMIRP  =  PERCENTAGE  FORM  OF  SMIR 

SNOW  =  INDICATOR  OF  SNOW  COVERAGE 

STDYR  =  STANDARD  YEAR  RADIATION 

STN  =  STATION  NUMBER 

STR  =  STATE  TAX  RATE 

STRP  =  PERCENTAGE  FORM  OF  STR 

TAC  =  TOTAL  COLLECTOR  AREA 

TAVC  =  PROPERTY  TAX  RATE  BASED  UPON  ASSESSED  VALUE  -  CONVENTIONAL 

TAVS  =  PROPERTY  TAX  RATE  BASED  UPON  ASSESSED  VALUE  -  SOLAR 

TBAR  =  EFFECTIVE  TAX  RATE 

TBTPR  =  TIME  BETWEEN  OUPUTS 

TCOST  =  TOTAL  COST  OF  AUXILIARY  HEAT  TO  DRYER  INLET  FOR  THE  STUDY 

TCR  =  TAX  CREDIT  RATE 

TCRP  =  PERCENTAGE  FORM  OF  TCR 

TDIC  =  DISCOUNTED  VALUE  OF  INCOME  TAX  DEDUCTIONS  ON  THE  INTEREST  - 

CONVENTIONAL 

TDIS  =  DISCOUNTED  VALUE  OF  INCOME  TAX  DEDUCTIONS  ON  THE  INTEREST  - 

SOLAR 

TDQAUX  =  DAILY  TOTAL  HEAT  PROVIDED  BY  AUX  HEATER  TO  DRYER  INLET 

THIN  =  INLET  OR  INITIAL  GRAIN  TEMPERATURE  (F) 

TIC  =  TEMPERATURE  INTO  COLLECTOR 

TIDD  =  DRY  BULB  TEMPERATURE  INTO  THE  DRYER 

TIDR  =  TEMPERATURE  INTO  DEHUMIDIFIER,  REGENERATION  SIDE 

TIDW  =  WET  BULB  TEMPERATURE  INTO  THE  DRYER 

TIME  =  TIME  (Ctf-WULATIVE)  (HR) 

TIN  =  INLET  AIR  TEMPERATURE  (F) 

TMIN  =  MINIMUM  TEMP  OF  DESICCANT  REGENERATION  AIR  (OPTION  *2> 

TOC  =  TEMPERATURE  OUT  OF  COLLECTOR 

TODP  =  OUTLET  TEMPERATURE  OF  DEHUMIDIFIER,  PROCESS  SIDE 

TODR  =  OUTLET  TEMPERATURE  OF  DEHUMIDIFIER,  REGENERATION  SIDE 

TOXP  =  OUTLET  TEMPERATURE  OF  SENSIBLE  HEAT  EXCHANGER,  PROCESS  SIDE 


TOXR 

TGAUX 

TQU 

TGDEHU 

TQDUCT 

TQHEX 

TQSURP 

TRD 

TT 

U 

VC 

VOID 
VS 

UAM 

WATER 

WEAT 

WDIR 

WIDP 

WIDR 

WODP 

UIODR 

WOXP 

WOXR 

WVEL 

XE 

XMAVE  ■ 

XME 

XMEND  i 

XMO 

YR 
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=  OUTLET  TEMPERATURE  OF  SENSIBLE  HEAT  EXCHANGER,  REGENERATION  SIDE 

=  TOTAL  AUXILIARY  HEAT  SUPPLIED,  KJ,  TO  DEHUMIDIFIER  REGENERATION  SIDE 

■  TOTAL  HEAT  GAIN  IN  COLLECTOR,  KJ 

=  HEAT  TRANSFERRED  IN  DEHUMIDIFIER,  KJ 

=  HEAT  LOST  FROM  DUCTING,  KJ 

=  HEAT  TRANSFERRED  IN  HEAT  EXCHANGER,  KJ 

=  TOTAL  SURPLUS  SOLAR  HEAT  AVAILABLE,  KJ 

=  INLET  TEMPERATURE  TO  AUXILIARY  HEATER,  REGENERATION  SIDE 

=  TOTAL  TIME  (IN  HOURS)  OF  STUDY 

=  EFFECTIVE  CONDUCTANCE  BETWEEN  DESICCANT  AND  AIR  STREAM 

=  RATIO  OF  ASSESSED  VALUATION  OF  THE  SYSTEM  IN  FIRST  YEAR  TO  THE  INITIAL 

IWESTMENT  OF  THE  SYSTEM  -  CONVENTIONAL 
'  VOID  FRACTION  OF  DESICCANT  MATERIAL 
=  RATIO  OF  ASSESSED  VALUATION  OF  THE  SYSTEM  IN  FIRST  YEAR  TO  THE  INITIAL 

INVESTMENT  OF  THE  SYSTEM  -  SOLAR 
:  AMBIENT  HUMDITY  RATIO 
:  AMOUNT  OF  H20  REMOVED  (CUMMULATIVE) 
'  WEATHER  CONDITIONS,  SNOW  ,  RAIN,  ETC. 
'  WIND  DIRECTION 

;  INLET  HUMIDITY  DEHUMIDIFIER,  PROCESS  SIDE 
;  INLET  HUMIDIPr  DEHUMIDIFIER,  REGENERATION  SIDE 

OUTLET  HUMIDITY  DEHUMIDIFIER,  PROCESS  SIDE 
'  OUTLET  HUMIDITY  DEHUMIDIFIER,  REGENERATION  SIDE 

OUTLET  HUMIDITY  HEAT  EXCHANGER,  PROCESS  SIDE 

OUTLET  HUMIDITY  HEAT  EXCHANGER,  REGENERATION  SIDE 

WIND  VELOCITY 

SENSIBLE  HEAT  EXCHANGER  EFFICIENCY 

AVERAGE  MOISTURE  CONTENT  IN  X-DIRECTION   (DECIMAL, DB> 

EQUILIBRIUM  MOISTURE  CONTENT   (DECIMAL, DB) 

FINAL  DESIRED  MOISTURE  CONTENT   (DECIMAL,  DB) 

INLET  OF  INITIAL  MOISTURE  CONTENT  (DECIMAL,  DB) 

YEAR 
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APPENDIX  I 

SAMPLE  OUTPUT 

OF 

DRYER  SIMULATION 
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ABSTRACT 

The  -feasibility  of  drying  grain  using  dehumidified  air  was  studied.  The 
dehumidified  air  was  produced  by  a  solar  energy  powered  sorption 
dehumidif ier.  A  grain  drying  model  was  coupled  with  an  existing  computer 
model  of  the  sorption  dehumidif ier .  Drying  grain  with  dehumidified  air  was 
compared  to  drying  grain  with  ambient  air,  solar  heated  air,  and  auxiliary 
heated  air.  The  air  that  was  dehumidified  was  found  to  dry  faster  than  air 
that  had  been  heated  only.  The  cost  per  bushel  was  found  to  be  about  tl.46 
for  the  conventional  system.  The  dehumidified  air  costs  were  *1.17/bu.  and 
the  dehumdif ied/auxil iary  heated  air  cost  »1.28/bu.  The  cost  per  bushel  of 
the  solar/auxiliary  heated  drying  air  was  tl.li.  A  study  of  the  life  cycle 
costs  of  the  solar  energy  powered  sorption  dehumdifier  showed  the  costs  of 
such  a  system  to  be  about  $2000.00  greater  than  a  conventional  drying 
system  that  uses  liquified  petroleum  gas.  A  grain  drying  system  that  uses 
a  solar  collector  to  provide  part  of  the  energy  needed,  was  found  to  have  a 
life  cycle  cost  of  about  *350  more  than  the  conventional  system. 


